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ABSTRACT 


A  lidar  was  designed  and  manufactured  at  the  Naval  Postgraduate  School,  Monterey, 
California,  to  provide  range  information  regarding  atmospheric  features,  such  as  clouds.  It 
is  further  planned  to  integrate  the  lidar  with  the  Naval  Postgraduate  School’s  Infrared 
Search  and  Target  Designation  system  (NPS-IRSTD)  at  some  future  date.  The  NPS- 
IRSTD  uses  two  vertical  infrared  linear  focal  plane  arrays  for  target  detection,  and  target 
direction  can  be  determined  very  accurately  but  the  system  does  not  provide  any  useful 
range  information.  The  lidar  was  proposed  as  the  solution  for  this  shortcoming.  The  lidar 
used  a  frequency-doubled  Nd:YAG  laser  which  had  an  energy  output  of  2  millijoules.  The 
laser  beam  was  expanded  to  17.75  inches  using  a  Dall-Kirkham  telescope  to  operate  within 
laser  safety  limitations.  The  theoretical  analysis  of  the  “Klett”  method  for  the  inversion  of 
lidar  returns  was  derived  and  a  MATLAB  program  was  written  to  demonstrate  the  process. 
A  daytime  and  nighttime  maximum  range  equation  for  the  lidar  was  developed.  The 
considerations  for  integrating  the  lidar  with  the  NPS-IRSTD  is  listed  and  a  solution  is 
proposed  to  obtain  the  mean  extinction  coefficient  along  the  path  in  the  infrared  spectrum 
using  the  lidar  inversion  extinction  coefficient  profile  at  532  nanometers. 
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I.  INTRODUCTION 


A.  BACKGROUND 

In  January  1984,  the  Naval  Postgraduate  School  (NPS)  received  the  Advanced 
Development  Model  of  the  AN/SAR-8  Infrared  Search  and  Target  Designation  (IRSTD) 
System  which  was  then  modified  to  bypass  the  background  normalization  circuitry  and 
designated  the  NPS-IRSTD.  The  NPS-IRSTD  is  a  360-degree  passive  surveillance  system 
for  the  detection,  tracking  and  identificadon  of  muluple  low-flying  targets.  The  ongoing 
modificadons  to  the  system  have  been  the  subject  of  many  theses  and  presentauons  in  the 
area  of  real  time  imaging  and  background  analysis  [Ref.  1],  The  NPS-IRSTD  uses  two 
linear  focal  plane  arrays  which  operate  in  the  3  to  5  micron  range  but  the  scene  image, 
presented  after  processing,  does  not  provide  any  useful  range  informadon  for  target 
declaration  to  a  combat  system. 

In  order  to  provide  the  needed  range  informadon  for  target  declaration,  one  of  two 
approaches  could  be  used,  passive  or  active.  One  passive  approach  would  use  a  three- 
color  method  and  would  require  three  separate  wavelength  regions  from  which  range  could 
be  determined  [Ref.  2].  The  NPS-IRSTD  would  have  to  be  modified  with  narrowband 
filters  to  provide  two  of  the  three  wavebands;  but,  to  provide  the  third,  an  additional 
detector  array  would  have  to  be  inrtnMed.  An  alternative  passive  approach  would  be  a  two¬ 
sensor  triangulation  but  this  would  require  an  additional  IRSTD  [Refs.  3  and  4],  The 
second  method  would  be  the  use  of  an  active  range-measurement  device,  such  as  a  radar  or 
laser.  The  laser,  used  in  this  capacity,  would  have  an  advantage  over  a  radar  because  it 
could  be  directed  to  irradiate  a  specific  area  of  interest  as  chosen  by  an  NPS-IRSTD 
operator  to  obtain  a  target  range.  A  disadvantage  of  the  laser,  besides  being  detectable  by 
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the  target,  is  that  it  is  not  an  all-weather  device  and  is  severely  attenuated  by  atmospheric 
losses  such  as  fog  and  lain.  However,  this  may  not  be  as  much  of  a  problem  as  first 
thought  since  the  NPS-IRSTD  operates  in  the  same  general  region  as  the  laser  and  would 
be  restricted  by  the  same  conditions. 

B.  RESEARCH  OBJECTIVES 

The  use  of  a  laser  as  a  rangefinder  is  not  a  new  concept  and  its  capabilities  extend 
into  the  radar  and  remote  sensing  fields.  It  would  be  useful  at  this  point  to  introduce  two 
new  terms  in  order  to  differentiate  between  their  uses.  The  term  “Laser  Radar”  implies  that 
the  laser  is  used  in  the  capacity  of  a  radar  to  detect,  track  and  image  targets.  The  term 
“LIDAR”,  which  is  an  acronym  for  Light  Detection  And  Ranging,  implies  that  the  laser  is 
used  for  atmospheric  remote  sensing  and  the  target  is  now  on  the  molecular  level.  In  both 
of  these  fields,  considerable  effort  has  been  expended  to  understand  the  concepts. 

The  lidar,  as  it  is  presently  planned  for  use  as  an  adjunct  to  the  NPS-IRSTD,  will  be 
required  to  provide  range  data  to  the  NPS-IRSTD  as  well  as  to  determine  the  atmospheric 
extinction  coefficient  along  the  path  to  the  target.  The  importance  of  this  extinction 
coefficient  information  would  be  realized  in  the  image  processing  of  the  NPS-IRSTD  image 
to  correct  for  atmospheric  propagation  losses  and  could  lead  to  faster  target  identification. 
Therefore,  the  primary  objective  of  this  thesis  was  to  manufacture  locally  a  lidar  which,  in 
the  future,  will  be  integrated  into  the  NPS-IRSTD.  The  range  and  remote  sensing 
capabilities  of  this  lidar  were  studied  and  documented  but  it  was  intended  to  restrict  the 
remote  sensing  capabilities  to  a  measurement  demonstration  of  the  atmospheric  extinction 
coefficient  at  this  time. 
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A  number  of  design  considerations  in  the  manufacture  of  the  lidar  were  taken  into 
account,  such  as  wavelength,  range  and  safety.  The  safety  consideration  was  one  of  the 
driving  concerns  due  to  die  op  .uonal  area  of  the  proposed  lidar,  which  will  be  on  top  of 
Spanagel  Hall  at  the  Naval  Postgraduate  School.  The  use  of  an  unsafe  laser  in  the 
surrounding  skv  is  not  permitted.  The  subsidiary  objective  of  this  thesis  was  to  detail  all 
the  desif  .1  considerations  and  develop  an  eye-sale  lidar  based  on  materials  available  at  the 
school. 

The  integration  of  the  lidar  into  the  NPS-1RSTD  system  was  not  an  objective  for  this 
thesis  but  a  preliminary  system  study  was  undertaken  as  a  second  subsidiary  objective. 

C.  NPS-IRSTD  SYSTEM  OVERVIEW 

The  NPS-IRSTD  system  consists  of  a  scanner  assembly,  a  twelve-channel  tape 
recorder  and  a  80386  computer  with  a  frame  grabber  data  acquisition  board  to  provide  real¬ 
time  scene  display.  It  has  a  field  of  view  of  360  degrees  by  10.5  degrees  and  operates  in 
the  3  to  5  micron  wavelength  region.  The  system  rotation  rate  is  0.5  Hertz  and  the  angular 
location  information  is  obtained  from  a  position  sensor  mounted  in  the  base  of  the  scanner. 
An  azimuthal  sample  width  of  0. 1047  milliradians  is  achieved  using  a  digitization  rate  of 
60,000  samples  per  revolution.  |  Ref.  1  ] 

The  scanner  assembly  consists  of  a  14  inch  F/l  Schmidt  telescope  with  a  10  inch 
entrance  aperture.  In  the  telescope's  focal  plane  ore  two  independent  vertical  linear  arrays 
of  indium  antimonide  (InSb)  detectors.  Each  linear  array  has  90  elements  and  each  element 
has  an  instantaneous  field  of  view  of  2.0  by  0.3  milliradians.  The  output  signals  from  each 
90-element  array  are  multiplexed  into  six  channels  and  are  then  converted  into  digital  form 
for  processing. 
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The  digital  data  is  fed  to  a  twelve-channel  tape  recorder  for  storage  at  a  data  rate  of  5 
megabytes  per  second  from  die  multiplexed  180  sensor  element  outputs.  Real-time  scene 
imaging  is  achieved  by  using  a  H0386  computer  with  a  frame  grabber  data  acquisition  board 
and  software  developed  locally.  |Ref.  1 1 

D.  LIDAR  TECHNOLOGY 

Lidars  are  used  in  the  study  of  the  atmosphere  and  the  ocean.  The  possible 
atmospheric  applications  for  lidar  are  the  measurements  of  extinction  profiles,  major  and 
minor  gas  densities,  temperature  and  humidity  profiles,  wind  speed,  cloud  composition, 
and  pollution.  The  oceanic  applications  include  measurements  of  surface  pollution,  sea 
surface  temperature  and  depth.  There  is  no  one  technique  alone  which  will  encompass  all 
applications  but  there  is,  instead,  a  range  of  application-dependent  approaches,  such  as 
molecular  fluorescence,  differential  absorption,  Raman  or  inelastic  scattering,  Mie/Rayleigh 
or  elastic  scattering  and  doppler  shift.  [Ref.  5] 

There  are  a  number  of  commercial  companies  which  manufacture  lidar  systems  but 
they  are  mainly  designed  to  meet  the  needs  of  the  customer’s  application.  Therefore,  it 
appears  that  the  lidar  is  still  a  specialty  tool  which  is  used  mainly  by  the  scientific 
community. 

E.  DESIGNED  LIDAR  SYSTEM  OVERVIEW 

The  lidar  designed  in  this  thesis  was  basic  and  limited  to  the  determination  of  the 
distance  to  clouds  based  on  measured  path  extinction  profiles.  Other  lidar  applications  are 
dependent  on  additional  technologies  which  were  not  available  at  the  time  of  manufacture 
but  may  be  considered  for  future  development. 
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The  designed  lidar  system  consists  c  a  transmitter,  a  receiver  and  electronic 
processing  equipment.  In  the  transmitter,  a  frequency-doubled  Nd:YAG  laser  with  a  pulse 
width  of  4  nanoseconds  and  variable  energy  output  was  used.  For  eye  safety,  the  laser 
beam  was  expanded  to  approximately  18  inches  in  diameter  before  propagation  in  the 
atmosphere.  The  laser  safety  aspects  are  detailed  in  Chapter  III.  The  lidar  receiver  consists 
of  a  lens  assembly  which  is  used  to  collect  the  backscattered  laser  pulse,  a  narrowband 
filter  and  a  photomultiplier  tube  detector.  The  electronic  processing  equipment  was  not 
fully  developed  for  this  thesis  but  the  output  of  the  photomultiplier  tube  was  observed  on  a 
digital  oscilloscope,  digitized  and  downloaded  to  a  computer  for  storage.  The  data  was 
later  analyzed  off-line  to  determine  the  atmospheric  extinction  coefficient. 

This  thesis  is  organized  as  follows: 


1 .  Chapter  II  is  an  introduction  to  lidar  theory  and  the  Klett  lidar  inversion  method. 

2.  Chapter  III  is  the  design  of  a  lidar  including  laser  safety  considerations  and 
equipment  specifications. 

3.  Chapter  IV  is  the  analysis  of  the  actual  data  obtained  from  the  designed  lidar. 

4.  Chapter  V  is  the  system  study  concerning  the  integration  of  the  lidar  into  the  NPS- 
IRSTD  system. 

5 .  Chapter  VI  is  a  discussion  of  the  conclusions  with  recommendations  for  future  lidar 
work. 

6.  Appendix  A  is  the  mathematical  analysis  upon  which  the  Klett  lidar  inversion  method 
is  based. 

7 .  Appendix  B  contains  theoretical  lidar  inversions  and  analyses  based  on  typical 
atmospheric  extinction  profiles  using  the  Klett  lidar  inversion  method. 

8.  Appendix  C  contains  the  detailed  experimental  data,  the  inversion  computer  program 
and  the  inverted  atmospheric  extinction  coefficient  profiles. 
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II.  LIDAR  THEORY 


A.  INTRODUCTION 

Over  the  years,  lidar  has  developed  into  a  useful  tool  for  remote  sensing  of  the 
atmosphere.  This  has  been  due  to  the  effort  of  many  dedicated  researchers  and  the 
invention  of  the  laser.  In  the  early  1930’s,  searchlights  were  first  employed  as  lidars  to 
assess  the  vertical  density  profile  of  the  upper  atmosphere  using  scattering  techniques  [Ref. 
6].  The  requirement  for  a  monochromatic  light  source  was  identified  when  the  scattering 
mechanisms  in  the  atmosphere  were  determined  to  be  wavelength  dependent.  With  the 
invention  of  the  laser,  researchers  began  to  realize  its  potential  for  remote  sensing  of  the 
atmosphere  by  using  previously  developed  searchlight  techniques.  Remote  sensing  of  the 
atmosphere  was  seen  to  hold  the  possibility  of  fast  and  economical  measurements 
compared  to  existing  methods  using  weather  balloons.  The  first  actual  demonstration  of  a 
working  lidar  appears  to  have  been  conducted  in  June  and  July  1963  by  Fiocco  and 
Smullin’s  use  of  a  ruby  laser  to  detect  backscatter  from  the  atmosphere  to  heights  of  140 
km  [Ref.  7]. 

In  its  simplest  form,  a  lidar  consists  of  an  optical  transmitter,  receiver  and  electronic 
processing  equipment,  A  typical  pulsed  optical  transmitter  requires  a  Q*switched  laser,  of 
suitable  energy  output  and  short  pulse  duration,  and  beam  expanding  optics.  The  receiver 
is  typically  composed  of  a  telescope  to  collect  the  backscattered  laser  pulse  and  a  sensor. 
The  electronic  processing  equipment  mathematically  transforms  or  inverts  the  time 
dependent  output  of  the  detector  to  obtain  the  atmospheric  data.  The  laser  and  optical 
receiver  assemblies  are  not  necessarily  mutually  exclusive;  there  are  lidars  in  service  in 
which  both  systems  share  the  same  telescope  housing. 
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B  .  LIDAR  EQUATION 


The  lidar  equation  is  the  comer  stone  to  the  use  of  the  lidar  and  can  be  derived  using 
the  Mie  theory  of  scattering  [Ref.  8|.  The  equation  defining  the  backscattered  power 
received  by  a  lidar  receiver  from  a  pulsed  laser,  assuming  single-scattering  and  omitting 
wavelength  dependence,  can  be  expressed  as: 


P(r)  = 


p(r)  exp 


CT(r')  dr’ 


(2-1) 


where 


P(r)  is  the  power  received  by  the  LIDAR  receiver  at  time  t, 

Pc  is  the  LIDAR  transmitted  power, 
c  is  the  velocity  of  light, 
t  is  the  duration  of  the  laser  pulse, 

Ao  is  the  effective  aperture  of  the  LIDAR  receiver, 

B(r)  is  the  atmospheric  volume  backscattering  coefficient  at  range  r, 
r  is  the  range  defined  by  c(  t-to)/2.  and 

a(r)  is  the  atmospheric  volume  extinction  coefficient  at  range  r. 


The  atmospheric  volume  backscatter  coefficient,  B(r),  is  the  backscattered  intensity 
per  unit  incident  intensity  per  unit  solid  angle  per  unit  path  length  (kilometers-1  • 
steradian-1)  and  characterizes  the  fraction  of  the  signal  scattered  back  along  the  propagation 
path  [Ref.  8],  B(r)  is  dependent  on  the  size  parameter  (which  determines  the  type  of 
radiative  scattering),  scattering  efficiency  and  the  number  density  of  the  scatterers  in  a  given 
volume.  The  size  parameter  is  defined  as  x  -  2  k  a  /X,  where  X  is  the  wavelength  of 
interest  and  a  is  the  radius  of  the  scattering  molecule  in  the  medium.  Therefore.  B(r)  is  a 
function  of  the  radiating  wavelength  and  size  of  molecular  scatterers  as  well  as  range.  A 
B(r)  range  dependence  is  introduced  because  the  number  density  of  molecular  scatterers  can 
vary  along  the  lidar  path  and  influence  the  magnitude  of  the  backscattered  laser  energy. 
B(r)  can  be  approximated  by  neglecting  the  diffuse  scattering  by  secondary  scatterers. 
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which  is  equivalent  to  assuming  that  the  backscattered  power  is  primarily  due  to  the  direct 
scattering  of  the  laser  radiation.  This  is  called  the  single  scattering  approximation;  it  may 
not  be  valid  in  very  dense  fog  [Ref.  9], 

The  atmospheric  volume  extinction  coefficient,  o(r),  defines  the  rate  at  which  a  signal 
attenuates  as  it  propagates  along  a  given  path,  expressed  in  kilometers1.  o(r)  is  dependent 
on  the  absorption  and  scattering  in  the  path  and  can  be  expressed  in  a  complete  form  as 
[Ref.  10]: 

o(r,X)  =  K(r,X)  +B(r,X,6)  (2-2) 


where 

v(r,X)  is  the  atmospheric  absorption  coefficient,  and 

B(rA,0)  is  the  atmospheric  scattering  coefficient  from  which  B(r,X)  is  derived. 

The  wavelength  dependence  is  clearly  stated  in  Equation  (2-2),  but  for  simplicity  this 
dependence  will  be  omitted  due  to  lidar  employment  of  near  monochromatic  laser  light. 
The  K(r)  term  can  be  further  broken  down  as  the  sum  of  the  aerosol  and  molecular 
absorption  coefficients,  while  B(r,0)  can  be  rewritten  as  the  sum  of  the  aerosol  and 
molecular  (Mie  and  Rayleigh)  scattering  coefficients.  Depending  on  the  spectrum  of 
interest,  there  are  regions  of  dominance  of  various  absorbers  and  scatterers.  In  the 
ultraviolet  spectrum,  the  large  ozone  molecular  K(r)  coefficient  dominates  the  value  of  o(r) 
while  in  the  visible  spectrum,  Baerosoi(r,0)  and  Bmoiecuiar(r,0)  strongly  influence  the  value 
for  o(r).  In  the  infrared  spectrum,  a(r)  is  dominated  by  the  various  molecular  absorption 
lines  of  which  water  and  carbon  dioxide  are  the  principal  absorbers;  this  directly  produces 
the  optical  windows  in  this  region.  The  range  dependence  is  given  because  the  number 
densities  of  the  scatterers  and  absorbers  can  vary  along  the  path  as  described  in  the 
discussion  on  B(r).  B(r,0)  can  be  related  to  B(r)  by  setting  the  scattering  observation  angle, 
0,  to  180  degrees  (backscatter). 
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It  has  been  observed  in  visible  and  near-infrared  experiments  that  a(r)  is  related  to 
U(r),  assuming  a  homogeneous  atmosphere,  by  the  empirical  power  law  relationship:  [Ref. 
11] 

P(r)  =  Constant  o(r)k  12-3) 

where 

k  is  range  independent  but  a  function  of  the  laser  wavelength  used  and  generally  in 
the  range  of  1.0  >  k  >  0.67,  and 

Constant  is  an  empirical  value  dependent  on  the  properties  of  the  scatterer 
irradiated. 

Therefore.  Equation  (2-3)  may  not  be  valid  in  those  areas  of  the  spectrum  where  K(r) 
dominates  or  strongly  influences  the  value  for  a(r),  such  as  in  the  8-12  micron  region. 
This  will  restrict  the  lidar’s  operational  spectrum  and  performance.  This  will  be  discussed 
in  the  last  section  of  this  chapter. 

C.  INVERSION  OF  THE  LIDA R  EQUATION 

The  lidar  equation  can  be  solved  for  the  volume  extinction  coefficient,  a(r),  by 
substituting  Equation  (2-3)  into  Equation  (2-1).  A  convenient  form  of  the  received  power, 
S(r)  =  In  (  r2  P(r)  ),  is  now  introduced  which  removes  the  range  dependence  from  the 
returned  signa.  and  expresses  the  result  in  a  logarithmic  form  for  analysis  purposes  [Ref. 
12].  The  lidar  equation  can  be  rewritten,  using  these  definitions  and  the  known  P(r0)  from 
a  reference  range  rQ  (the  minimum  range  for  the  overlap  of  the  receiver’s  and  transmitter’s 
fields  of  view): 


S(r)  -  S(r0)  =  In 


a(r’)  dr'  . 


(2-4) 
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The  lidar  equation  can  be  further  reduced  by  differentiating  both  sides  with  respect  to 


range.  The  resulting  differential  equation  is  as  follows: 


dS(r) 

dr 


dft(r) 

dr 


-  2a(r)  . 


(2-5) 


The  expression  for  B(r),  Equation  (2-3),  can  be  substituted  into  Equation  (2-5)  and 
results  in  the  following  differential  equation: 


dS(r)  _  k  dq(r) 
dr  o(r)  dr 


-  2o(r)  . 


(2-6) 


The  differential  Equation  (2-6)  is  nonlinear  because  of  the  kJa  da/dr  term  but  has  the 
form  of  a  Bernoulli  differential  equation.  It  can  be  shown,  by  using  a  linearization 
transformation  (n  =  l/a(r)),  that  the  resulting  solution  to  the  differential  equation  is: 


o(r)  = 


_ exp[(S(r)  -  S(r0))/k] _ 

-7—  -  £  I  exp[(S(r’)-S(r0)yic]  dr¬ 
oll  0)  k  I 


(2-7) 


A  number  of  inversion  techniques  have  been  used  to  solve  Equation  (2-7)  [Refs.  11, 
13,  14,  15,  16,  17,  and  18]  and  it  is  not  the  purpose  of  this  section  to  describe  all  of  them. 
One  technique,  however,  has  become  the  standard  against  which  most  other  techniques  are 
compared;  it  is  the  “Klett  method"  [Ref.  1 1], 
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1 .  Klett  Lidar  Inversion  Method 

In  the  Klett  method,  instead  of  integrating  Equation  (2-7)  forward  to  solve  for 
the  extinction  coefficient  from  r  to  r„,  a  boundary  value  (<r(r,n))  for  the  extinction 
coefficient  is  determined  at  the  lidar 's  maximum  range  (rm).  Then  the  extinction  coefficient 
can  be  obtained  by  integrating  backwards  from  rm  to  r.  This  method  has  a  number  of 
advantages  such  as  the  received  signal  strength  being  attenuated  only  to  a  known  value 
(P(rm))  rather  than  starting  at  P(r„)  and  decaying.  With  sensor  systems,  there  will  always 
be  a  finite  limit  at  which  a  signal  can  be  detected  in  noise  and  the  Klett  method  is  better 
suited  to  a  practical  application.  The  resulting  Klett  equation  is  as  follows: 


o(r)  = 


_ exp  [(S(r)  -  S(rm))/k] _ 

ftm 

—1—  +  J  I  exp[(S(r')  -  S(rm)yk]  dr' 

<r(rm)  k  / 


The  determination  of  the  extinction  coefficient  boundary  value,  afim),  becomes 
very  important.  The  value  for  ct( rm)  can  be  approximated  by  relating  it  to  the  signal 
strength  S(rm).  Appendix  A  contains  the  derivation  of  this  relationship.  In  summary,  the 
relationship  can  be  obtained  through  a  variable  substitution  into  Equation  (2-7),  assuming 
that  the  average  extinction  coefficient  over  the  range  0  to  rm  is  equal  to  that  for  r0  to  rm 
[Ref.  19]  and  then  substituting  the  results  into  the  defining  expression  for  S(r).  The 
relationship  between  c( rin)  and  S(rm)  can  be  expressed  as: 


S(rm)  -  C 


+  ln(2(rm  -r°))  =  ln  (Qm)  -  - — 5a»—  In  ( 1  +  I  Q.m)  . 

'  K  '  \rm  *  r0J 
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Here  the  terms  C,  I  and  £2m  are  introduced  and  defined  as  follows: 


1 .  C  is  a  grouping  of  lidar  constants  from  the  expression  for  Slr^  defined  as: 

C  =  In  |p„  ^  A  Constant  J  (2-10) 

2 .  I  is  a  convenient  dimensionless  integral  term  defined  as: 


1  ~  (fin  *  fi 


f  rm 

o)'1  0 

Jto 


exp[(S(r’)  -  S(rm))/k]  dr' 


(2-11) 


3.  £2m  is  also  a  convenient  dimensionless  grouping  which  is  related  to  the  expression 
for  I.  £2m  can  also  be  used  as  a  measure  of  optical  depth  over  the  range  (r0,rm)  if  the 
average  extinction  along  the  path  is  equal  to  a(rm).  £2m  is  defined  as: 


ri  _  2  o(rm)  (rin  -  r0) 
m  ~  k 


(2-12) 


A  further  formal  simplification  results  from  setting  the  left  side  of  Equation  (2- 
9)  equal  to  Gm,  which  is  a  constant  term.  The  resulting  simplified  relationship  is  then: 


Gm  =  ln(£2m)---&_  In ( 1  +  I £2m)  . 

(fin  *  roJ 


(2-13) 


Since  Gm  is  a  known  quantity,  the  determination  of  £2m  becomes  simply  an 
iterative  root  finding  process.  In  Figure  1  is  the  graphical  representation  of  the  solution 
curve  for  Equation  (2-13)  where  the  y  axis  is  defined  as  In  (£2)  -  fin/(fin  -  r0)  In  (1  +  I  £2). 
When  the  solution  curve  intersects  with  the  value  for  Gm,  a  value  for  o(rm)  can  easily  be 
found  using  Equation  (2-12).  The  selection  of  the  root  for  £2m  is  ambiguous  unless  the 
range  and  visibility  are  considered.  The  solution  for  the  correct  root  for  £2m  is  determined 
by  looking  at  the  limiting  cases  of  high  and  low  visibility  conditions  and,  failing  these 
cases,  applying  a  default  algorithm.  A  summary  of  the  procedure  is  as  follows:  [Ref.  20] 
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1 .  First  try  the  High  Visibility  algorithm.  Using  the  assumption  that  the  extinction 
coefficient  over  the  range  0  to  r0  is  constant,  then  the  value  for  the  integral  in 
Equation  (2-1)  in  this  region  can  be  written  as  2  <j(r0)  rQ.  A  new  expression  for  Gm 
can  be  obtained  by  modifying  Equauon  (2-13)  for  the  total  path  0  to  rm. 


Gm  =  Gm  +  2  r„  a(r())  /k  =  In  (Qm)  -  In  ( 1  +  1  Qm)  . 
Equauon  (2-14)  can  be  rewritten  in  terms  of  with  the  resultant  form: 

=  [exp  (Gm  - 1)]  . 


(2-14) 


(2-15) 


The  value  of  is  obtained  iteratively  and  its  value  is  assumed  correct  if  exp(G'm)  > 
(I  +  0.01),  if  the  resulting  value  for  a(rm)  >  0.01  and  if  a(r0)/<J(rm)  <  50. 

2.  If  the  High  Visibility  algorithm  fails,  then  try  the  Low  Visibility  algorithm.  When 
operating  in  low  visibility  conditions,  it  can  be  shown  that  the  mean  value  of  the 
extincuon  coefficient  will  depend  primarily  on  the  magnitude  of  the  backscattered 
signal  which  influences  the  value  for  I.  A  simple  approach  is  used  in  which  the 
mean  of  the  extincuon  coefficient  is  set  equal  to  the  boundary  value,  afr^)  (rm  -  rQ). 
The  resulting  equauon  can  be  solved  iterauvely  for  Qm.  and  Cfim)  determined. 

Qm  =  ln(l  +IftJ  .  (2‘16) 


The  algorithm  gives  a  solution  if  I  >  1. 

3 .  If  the  previous  two  algorithms  fail,  then  apply  the  Default  algorithm.  The  default  is 
set  to  the  maximum  value  for  when  there  is  no  soluuon  to  the  Klett  equauon. 

o  _(rm-r„)  (2-17) 

**max  —  r  • 

r0  I 
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A  theoretical  analysis  of  lidar  returns  is  conducted  in  Appendix  B  using  typical 
extinction  coefficient  proxies.  The  effects  of  the  lidar  range  resolution  and  atmospheric 
noise  are  considered  in  the  analysis  for  their  effect  on  the  inversion  process.  The  lidar’s 
range  resolution  is  strictly  a  function  of  die  data  sampling  rate  used  to  digitize  the  lidar 
receiver’s  analogue  signal.  The  MATLAB  high-performance  numeric  computation 
software  is  used  for  the  analysis  which  produces  a  S(r)  file  from  a  specified  path  extinction 
coefficient  profile  and  inverts  the  S(r)  file  based  on  the  algorithms  outlined  above.  In 
some  cases,  the  Low  Visibility  algorithm  was  observed  to  be  unreliable,  which  was  also 
noted  by  KletL  The  simple  assumption  that  the  mean  of  the  extinction  coefficient  is  equal 
to  the  boundary  value.  Equation  (2-16),  is  invalid  when  there  is  a  large  dynamic  range  of 
the  extinction  coefficient.  The  solution  for  low  visibility  cases  occurs  when  the  Gm  line 
and  Qm  curve  are  almost  parallel  and  any  analysis  errors  are  magnified  in  the  solution. 

D.  RECEIVER  THEORY 

The  operation  and  performance  of  the  lidar  receiver  are  as  important  as  the  lidar 
transmitter  itself.  In  this  section,  the  theoretical  performance  of  a  photomultipler  tube 
receiver  will  be  studied  and  characterized.  This  will  lay  the  ground  work  for  the  next 
chapter  concerning  system  design  . 

The  photomultiplier  tube  is  a  vacuum  tube  photoemissive  detector  which  consists  of  a 
cathode,  a  number  of  dynodes  and  an  anode  as  shown  in  Figure  2.  The  cathode  is  a 
photocathode  which  generates  free  electrons  when  struck  by  photons  of  sufficient  energy  to 
overcome  the  cathode  work  function.  The  generated  electrons  are  accelerated  towards  the 
anode  by  the  electric  field  between  the  anode  and  cathode  but  are  deflected  to  strike  a 
dynode.  The  dynodes  are  electrodes  at  successively  greater  potentials  and  are  made  up  of  a 
material  which  has  a  low  work  function.  The  free  electron  collision  with  the  dynode  results 
in  the  production  of  several  secondary  electrons.  These  electrons  are  then  directed  to  the 
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next  dynode  and  gains  are  achieved  in  the  order  of  106  through  successive  dynodes  and 
secondarv  •'  ctron  pK  -lion.  This  gain  can  be  characterized  by  G  =  5N,  where  G  is  the 
gain,  N  is  the  number  of  dynodes  and  5  is  the  number  of  secondary  electrons  per  impact 
which  is  a  function  of  interdynode  voltage  and  dynode  material. 

The  operational  use  of  photomultiplier  tubes  has  been  restricted  due  to: 

1 .  The  fragile  nature  of  vacuum  devices. 

2.  The  high  voltage  requirements  to  achieve  sufficient  gains  (100s  of  volts  per  stage). 

3 .  The  wavelength  restrictions  because  photocathode  spectral  response  is  limited  to 
visible  and  ultraviolet  wavelengths. 

4.  The  low  quantum  efficiencies,  typically  10-25%. 


Vacuum  Hum 


Figure  2.  Photomultiplier  tube  typical  layout  [Ref.  21]. 

The  photomultiplier  tube  has  continued  to  be  manufactured  over  the  years  to  meet  the 
needs  of  the  scientific  community.  Its  desirable  features  are: 
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1 .  The  high  internal  gains  which  allow  for  the  detection  of  very  weak  signals. 

2 .  The  low  internal  noise. 

3.  The  last  response  times.  The  typical  response  times  for  photomultiplier  tubes  are  in 
the  order  of  2  nanoseconds. 


The  sources  of  noise  in  the  photomultiplier  tube  are: 

1 .  Shot  noise  (ifsj).  This  is  noise  due  the  statistical  generation  process  of  electrons  by 
the  cathode.  The  rms  value  of  the  shot  noise  current  can  be  written  as  [Ref.  22): 


iN  =  V 2q  (is  +  id  +  »b  )  Af  G2  (2-18) 

where 

is  the  average  signal  cathode  current  due  to  the  lidar  backscattered  photons, 
id  is  the  dark  current  (the  thermally  generated  current  when  no  photon  flux  is 
present  at  the  cathode), 

ib  is  the  average  cathode  current  due  to  the  background  radiance 
(discussed  in  the  following  paragraph), 

Af  is  the  electrical  system  bandwidth,  and 
G  is  the  photomultiplier  tube  gain. 

2.  Background  noise  (ib).  This  is  the  shot  noise  component  due  the  background  scene 
which  can  be  limited  by  the  choice  of  the  receiver  field  of  view  and  narrowband 
optical  filters.  A  discussion  of  this  noise  source  component  is  felt  to  be  important 
because  it  will  be  the  limiting  factor  for  daylight  photomultiplier  operation.  The  total 
generated  cathode  current  is  due  to  the  summation  of  the  signal  and  average 
background  photon  flux.  These  currents  are  not  necessarily  independent  because  the 
transmitted  laser  pulse  will  excite  the  atomic  molecules  along  the  path.  The  molecular 
excited  state  population  density  is  increased  with  each  laser  pulse  and  there  is  a 
transition  time  required  to  return  to  the  original  state.  The  transition  process  can  be  a 
radiative  one,  thereby  increasing  the  total  background  photon  flux.  In  this  analysis, 
this  effect  will  be  ignored  because  it  mainly  applies  to  high  energy  lasers.  A  general 
expression  for  the  average  background  cathode  current,  assuming  no  increased 
background  flux,  can  be  written  as:  [Ref.  22] 


ib 


Xqri(X)  Eb(K) 

h  c  Td 


(2-19) 
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where 

X  is  the  wavelength  of  background  radiation. 
r|(X)  is  the  photomultiplier  tube  quantum  efficiency, 

Eb(X)  is  the  background  spectral  energy, 
c  is  the  speed  of  light, 
h  is  Planck’s  constant,  and 

tj  is  the  detection  interval  which  is  different  from  the  laser  pulse  duration. 


The  background  energy  can  be  written  in  terms  of  the  receiver’s  narrowband  filter 
and  solid  angle  assuming  that  the  individual  terms  which  influence  the  background 
energy  magnitude  are  constant  over  the  wavelength  band  interval: 


Eb(X)  =  AX  Sb(X)  t„(X)  flo  Ao  Td 


(2-20) 


where 

AX  is  the  wavelength  interval  of  the  narrowband  filter. 

Sb(X)  is  the  background  spectral  radiance  (This  term  can  be  derived 
from  the  radiative  transfer  equation  in  the  visible  spectrum  where  the 
principal  scattering  mechanism  is  Rayleigh  scattering), 
tQ(X)  is  the  receiver  optics  transmittance, 

£2o  is  the  receiver  solid  angle, 

Ao  is  the  receiver’s  effective  aperture,  and 
xd  is  the  detection  interval. 

The  resultant  average  background  cathode  current  can  be  written  as  follows: 


X  q  q(X)  AX  Sb(X)  xp(X)  fl0  A» 
h  c 


(2-21) 


3.  Johnson  or  thermal  noise  (ij).  This  is  the  noise  generated  by  the  thermal  excitation  of 
the  charge  carriers  in  a  resistor.  The  rms  value  of  the  Johnson  noise  current  can  be 
written  in  the  form: 

/4kTAf 

*j  rms  -  "y  ^  (2-22) 

where  R  is  the  equivalent  resistance  of  the  load  and  photomultiplier  tube  resistances. 
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4.  Dynode  noise.  This  is  the  noise  due  to  the  statistical  fluctuations  in  the  generation  of 
secondary  electrons  at  the  dynodes.  It  can  be  shown  that  the  noise  contribution  of  the 
first  dynode  dominates  the  noise  process  of  the  complete  dynode  amplification  but 
there  is  an  inverse  relationship  between  the  gain  of  the  first  dynode  and  the  signal-to- 
noise-ratio.  The  dynode  noise  term  can  be  neglected  assuming  that  the  gain  of  the 
first  dynode  is  high  enough  to  make  the  photomultiplier  tube  gain  a  noiseless 
process.  [Ref.  21] 


The  general  expression  for  the  power  signal-to-noise  ratio  (SNR)  across  a  load 
resistor  can  be  written,  using  Equations  (2-19)  and  (2-22),  as  follows  [Ref.  21]: 


SNR 


_ RLoad  *s _ 

R|,,ad(2q(i;+  id  +  i;)AfG2  +  i^p^) 


(2-23) 


A  number  of  limiting  conditions  can  be  applied  to  Equation  (2-23).  If  R  is  large 
enough  to  make  the  Johnson  noise  very  small  compared  to  the  shot  noise,  then  the  shot 
noise  will  dominate  the  noise  power  term  and  the  Johnson  noise  may  be  neglected.  The 
remaining  two  conditions,  daylight  and  nighttime  operation,  will  be  discussed  in  detail 
because  of  their  importance  in  determining  the  maximum  range  for  the  lidar  system  which 
will  become  apparent  later  in  the  next  two  subsections. 

1 .  Daytime  Lidar  Receiver  Operation 

If  the  photomultiplier  tube  is  used  for  daylight  operation,  then  the  average 
background  current  will  dominate  the  dark  current  and  signal  current  terms  for  the  shot 
noise  term.  Therefore,  the  dark  current  and  signal  current  terms  can  be  neglected  in  the 
analysis.  The  resultant  daylight  power  SNR  from  Equation  (2-23)  will  be: 


SNR 


is  R-Load 

2  q  ib  Af  G2  RLoad 


(2-24) 
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The  minimum  detectable  optical  power  is  defined  when  the  signal  power  across  the 
load  resistor  is  equal  to  the  noise  power,  or  SNR  =  1.  This  value  for  the  SNR  is  arbitrary 
but  represents  the  best  case  scenario  for  the  receiver  and  also  represents  the  Noise 
Equivalent  Power.  An  expression  for  the  minimum  detectable  incident  power  can  be 
determined  by: 

1 .  Noting  that  the  signal  current  at  the  cathode  can  be  expressed  in  terms  of  the  incident 
power  P(r)  which  is  range  (time)  dependent  due  to  scattering: 

X  q  r|(A.)  P(r)  (2-25) 

's  he 

2.  Substituting  Equations  (2-21 )  and  (2-25)  into  the  expression  for  the  SNR,  liquation 
(1-24). 

3.  Solving  the  resultant  expression  for  P(r),  with  SNR  =  1, 


P(r)  = 


I  2  h  c  AX  Sh(X)  T0(X)  fl0  Ap  Af 
V  \  r\(A.) 


(2-26) 


The  maximum  range  for  ihe  lidar  system  can  be  approximated  by  substituting  the 
expression  for  minimum  P(r)  (Equation  (2-26))  into  the  general  lidar  expression  (Equation 
(2-1))  and  solving  for  the  range,  r.  The  resulting  expression  for  the  maximum  daytime 
range  of  a  photomultiplier  lidar  system  is: 


(2-27) 


where 
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rma*  is  the  maximum  range  for  the  lidar  system  (meters), 

P0  is  the  lidar  system  transmitted  power  (watts), 

T  is  the  laser  pulse  duration  (seconds), 

li(r)  is  the  atmospheric  volume  backscattering  coefficient  at  range  rmax 
(meters- 1  •  steradians1), 
c  is  the  speed  of  light  (meters/sccond), 

Ao  is  the  receiver’s  effective  aperture  (meter2), 

X  is  the  wavelength  of  background  radiation  which  can  be  approximated  as  the  laser 
wavelength  (meters)  because  of  the  narrowband  optical  filter, 

T|(X.)  is  the  photomultiplier  tube  quantum  efficiency, 

CT(r)  is  the  atmospheric  volume  extinction  coefficient  (meters1), 
h  is  Planck’s  constant  (watts  •  seconds2), 

AX  is  the  wavelength  interval  of  the  narrow  band  filter  (microns), 

Sb(A,)  is  the  background  spectral  radiance  (watts  •  meter-2  .  steradians1  •  micron-1), 
T0(X)  is  the  receiver  optics  transmittance, 
is  the  receiver  solid  angle  (steradians),  and 
Af  is  the  electrical  system  bandwidth  (Hertz). 

2 .  Nighttime  Lidar  Receiver  Operation 

For  nighttime  operation,  the  signal  current  (is)  will  dominate  the  average 
background  current,  which  can  be  neglected.  The  resultant  nighttime  power  SNR  from 
Equation  (2-23)  will  be: 


SNR 


is  G2  ^Load 
2  q  is  At  G2  RLoad 


(2-28) 


or 


SNR  =  - i -  .  (2-29) 

2  q  is  Af 

A  simplification  of  Equation  (2-29)  is  possible  if  the  response  time  of  the 
photomultiplier  tube  is  faster  than  the  laser  pulse  width  used  in  the  lidar  transmitter.  Then 
is  over  the  laser  pulse  width  is  approximately  equal  to  is.  Then  Equation  (2-29)  can  be 
written  using  Equation  (2-25)  for  the  expression  for  is  as: 
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(2-30) 


SNR 


X  q  r)(A.)  P(r) 
2  q  h  c  Af 


An  expression  for  the  P(r)  can  he  obtained  by  setting  the  minimum  detectable 
optical  power  across  the  load  resistor  equal  to  the  noise  power  or  SNR  =  1 : 


P(r)  = 


2  h  c  Af 
Xti(X) 


(2-31) 


The  maximum  range  for  the  lidar  system  can  be  approximated  by  substituting 
the  expression  for  the  minimum  P(r)  (Equation  (2-31))  into  the  general  lidar  expression 
(Equation  (2-1))  and  solving  for  the  range.  The  receiver  optics  transmittance  (x0(A.))  must 
also  be  taken  into  account  when  developing  an  expression  for  P(r)  and  the  new  expression 
will  be  simply  Equation  (2-31)  multiplied  by  x0(^-).  The  resulting  expression  for  the 
maximum  nighttime  range  for  a  photomultiplier  lidar  system  is: 


! 

rrmax 

P„  x  Ao  x0(X)  X  ti(X)  P(r)  exp 

-2  1  a(r’)  dr' 

Jo 

(2-32) 


4  h  Af 


where 


fmax  is  the  maximum  range  for  the  lidar  system  (meters), 

P0  is  the  lidar  system  transmitted  power  (watts), 
x  is  the  laser  pulse  duration  (seconds), 

B(r)  is  the  atmospheric  volume  backscattering  coefficient  at  range  rmax 
(meters- 1  •  steradians-i), 
c  is  the  speed  of  light  (meters/second), 

Ao  is  the  receiver’s  effective  aperture  (meter2), 

X  is  the  wavelength  of  background  radiation  which  can  be  approximated  as  the  laser 
wavelength  (meters)  because  of  the  narrowband  optical  filter, 
t|(X)  is  the  photomultiplier  tube  quantum  efficiency, 
a(r)  is  the  atmospheric  volume  extinction  coefficient  (meters- •), 
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h  is  Planck’s  constant  (watts  •  seconds2), 
x0(X)  is  the  receiver  optics  transmittance,and 
Af  is  the  electrical  system  bandwidth  (Hertz). 

The  purpose  of  the  derivation  of  Equations  (2-27)  and  (2-32)  is  to  provide  the 
theoretical  basis  for  the  maximum  range  for  a  lidar  system  for  day  and  nighttime  operation. 
In  general,  the  maximum  nighttime  range  will  greater  than  the  daylight  range  because  of  the 
background  spectral  radiance  term  (St,$.))  present  in  Equation  (2-27)  which  increases  the 
optical  noise  in  the  lidar  system  and  decreases  the  range. 

E.  LIDAR  WAVELENGTH  SELECTION 

The  selection  of  the  lidar  wavelength  is  an  important  consideration  for  the  design  of 
any  lidar  system.  In  Section  B,  the  atmospheric  backscatter  coefficient  was  related  to  the 
atmospheric  volume  extinction  coefficient  in  order  to  solve  the  lidar  equation;  but  this 
relationship  was  wavelength  dependent.  The  main  options  for  an  affordable  lidar  laser 
wavelength  are  in  the  visible  (400-700  nanometers),  near-infrared  (700  nanometers-2.5 
microns)  and  far-infrared  (8-12  microns)  regions. 

In  the  visible  and  near-infrared  regions,  the  relationship  between  the  atmospheric 
backscatter  coefficient  and  the  atmospheric  volume  extinction  coefficient  has  been  shown  to 
be  constant  over  a  wide  range  of  size  parameters  [Ref.  23].  The  laser  beam  divergence  for 
lasers  in  these  regions  is  small  based  on  the  Gaussian  far-field  full-angle  beam  divergence 
formula,  1.27  times  laser  wavelength  divided  by  twice  the  beam  waist  diameter.  In  theory, 
a  larger  photon  flux  density  is  possible  due  to  the  smaller  beam  divergence  which  results  in 
a  larger  backscattered  signal  when  not  taking  into  account  wavelength  dependent  scattering 
mechanisms.  Therefore,  the  maximum  range  of  a  lidar  can  be  increased  by  using  a  laser 
with  a  smaller  beam  divergence.  When  considering  visible  and  near-infrared  lasers  for 
lidar  applications,  visible  lasers  will  have  an  advantage  based  on  beam  divergence  when 
the  laser  output  energies  are  the  same. 
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There  are  a  number  of  disadvantages  with  visible  and  near-infrared  lasers.  The  main 
consideration  is  eye  safety.  Lasers  in  the  visible  and  near-infrared  regions  can  be 
dangerous  to  the  retina  of  the  eve  and  energy  levels  must  be  kept  low  for  external  use  in 
urban  areas.  This  restricts  the  maximum  remote  sensing  range  of  the  lidar.  There  is  also 
the  problem  of  solar  background  radiance  which  can  wash  out  weak  lidar  signals  unless 
special  measures,  such  as  the  use  of  narrowband  filters,  are  taken  to  reduce  the  influence  of 
the  background  in  the  detector. 

In  the  far-infrared  region,  there  are  many  advantages.  The  primary  advantage  is  the 
issue  of  eye  safety  for  operation  in  urban  areas.  Those  lasers  operating  in  the  far  infrared 
region  are  live  orders  of  magnitude  safer,  based  on  Maximum  Permissible  Exposure  [Ref. 
24],  than  visible  lasers.  There  are  a  number  of  commercial  high  power  lasers  in  this  region 
which  are  readily  available,  such  as  the  carbon  dioxide  laser  (10.6  microns).  The  carbon 
dioxide  laser  has  one  of  the  best  efficiencies  and  highest  power  outputs  of  commercially 
available  lasers.  The  solar  radiance  problem  is  reduced  to  the  thermal  emission  of  the 
background  which  is  about  one  order  of  magnitude  smaller  than  the  visible  region 
background  radiance  [Ref.  22], 

As  with  visible  and  near-infrared  lasers,  the  far-infrared  laser  has  a  number  of 
disadvantages.  In  diis  region,  the  relationship  between  the  atmospheric  volume  backscatter 
and  atmospheric  volume  extinction  coefficient  is  not  as  linear  as  in  the  visible  region.  The 
magnitude  of  the  backscatter  coefficient  in  the  far  infrared  is  almost  two  orders  of 
magnitude  smaller  than  a  corresponding  visible  region  value  [Ref.  25],  The  result  for  far- 
infrared  lasers  is  a  one  hundred  fold  increase  in  the  pulsed  energy  output  required  to 
achieve  the  same  lidar  performance  as  a  visible  laser  (using  Equation  (2-1)),  but 
commercially  available  carbon  dioxide  lasers  can  provide  the  needed  energy  levels.  The 
final  disadvantage  of  lidar  operation  in  this  region  is  beam  divergence  which  is 
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approximately  1(X)  times  worse  than  in  the  visible  region  and  will  limit  the  maximum 
operating  range. 

The  choice  of  a  lidar  system  wavelength  is  not  an  easy  one  and  is  dependent  on  the 
employment  and  measurements  to  he  made.  In  this  thesis,  a  frequency-doubled  Nd:YAG 
laser  was  used  (532  nanometers)  mainly  due  to  availability  as  well  as  its  many  advantages 
over  far-infrared  lasers  as  discussed  above. 

The  next  chapter  will  deal  with  the  actual  design,  manufacture  and  calibration  of  a 
frequency-doubled  Nd:YAG  lidar  system. 
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III.  LIDAR  SYSTEM  DESIGN 


A.  OVERALL  DESIGN 

The  design  of  the  lidar  system  was  carried  out  for  the  maximum  possible  range 
subject  to  the  constraints  of  laser  availability  and  eye  safety.  The  designed  system  consists 
of  a  frequency-doubled  Nd:YAG  medical  laser,  a  17.75  inch  telescope  for  the  transmitting 
optics,  a  photomultiplier  tube  with  associated  optics  used  as  the  detector,  a  digital 
oscilloscope  for  data  collection  and  display  and  a  computer  for  data  acquisition,  storage  and 
off-line  processing.  A  schemauc  of  the  overall  lidar  design  is  featured  in  Figure  5.  The 
main  constraint  for  the  design  of  the  lidar  was  laser  safety  since  the  lidar  will  be  used  for 
atmospheric  studies  and  will  be  located  on  the  school  grounds. 

The  initial  lidar  design  was  based  on  a  large  area  Fresnel  lens,  but  in  testing  it  was 
found  to  be  unsuitable  for  beam  expansion.  The  Fresnel  lens  used  was  an  Edmund 
Scientific  experimental  Fresnel  lens  D34.747  which  was  approximately  28  inches  by  37 
inches  and  had  a  focal  length  of  54  inches  [Ref.  26].  The  40  lines  per  inch  on  the  lens  was 
not  sufficiently  high  enough  to  produce  a  good  quality  collimated  beam.  This  lidar  design 
could  not  progress  further  until  an  alternative  beam  expander  was  found.  A  search  of  the 
NPS  Physics  Department  revealed  an  old  telescope  which  was  used  in  the  final  design. 

B ,  LIDAR  TRANSMITTER 

1 .  Laser 

The  laser  used  in  the  lidar  transmitter  was  a  Nd:YAG  laser  model  MK  100 
manufactured  by  Kigre  Incorporated.  The  laser  was  later  modified  by  Kigre  Incorporated 
with  an  extra-cavity  frequency-doubling  crystal  to  provide  a  532  nanometer  wavelength 
output.  The  system  consisted  of  a  laser  head  and  separate  power  supply.  The  power 
supply  incorporates  a  control  unit  which  allows  the  selection  of  the  number  of  pulses  per 
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firing  (1-5),  the  energy  per  pulse  and  a  laser  energy  monitor  display.  A  summary  of  the 
laser’s  specifications  is  in  Tabic  I . 

After  the  frequency-doubling  process,  both  the  original  1.06  micron  and  the 
frequency-doubled  532  nanometer  wavelengths  were  present  in  the  laser  output.  As 
discussed  in  Chapter  2,  the  visible  wavelength  was  desired  due  the  strong  relationship 
between  o(r)  and  (3(r).  Therefore,  it  was  necessary  to  remove  the  1.06  micron  beam  from 
the  main  laser  beam.  It  was  noted  that  the  laser’s  energy  density  was  not  the  classic 
Gaussian  profile  as  expected  and  there  was  an  off-axis  laser  artifact  present  with  enough 
energy  to  produce  a  bum  mark  in  carbon  paper,  as  shown  in  Figure  6.  The  discovery  of 
the  artifact  was  important  as  it  required  removal  for  laser  safety  reasons.  The  Nd:YAG 
laser’s  external  optics  were  inspected  by  the  author  and  nothing  was  found  which  could 
suggest  the  cause  of  the  laser  beam  energy  profile.  It  was  therefore  assumed  to  be  a 
possible  misalignment  of  the  frequency-doubling  crystal  which  was  sealed.  The  only 
option  to  correct  this  problem  was  to  return  the  laser  to  Kigre  Incorporated  for  repair,  but 
this  would  have  imposed  unacceptable  delays  in  the  completion  of  the  thesis.  It  was 
therefore  decided  to  complete  the  thesis  and  investigate  the  possibilities  for  the  Nd:YAG 
laser  repair  at  a  later  date. 

It  was  observed,  after  a  large  number  of  laser  firings,  that  the  energy  output  of  the 
Nd:YAG  laser  varied  from  shot  to  shot  from  those  measured  in  Table  1.  This  should  not 
present  a  problem,  since  only  the  maximum  energy  values  will  be  used  in  any  laser  safety 
calculations. 
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TABLE  1.  Nd:YAG  LASER  SPECIFICATIONS  [Ref.  27] 


Characteristic 

Value 

Comment 

Maximum  Energy  Output 
(5  pulses  per  shot) 

mean  =  9.126  millijoules 
a  =  1 .587  millijoules 

Provided  by  Kigre 
Incorporated  for  a  532 
nanometer  100  snot  test. 

Beam  Diameter 

4.0  ±  0.2  millimeters 

Determined  by 
experimentation  using  laser 
bum  marks. 

Half  Angle  Beam 

Divergence 

0. 165  ±  0.001  milliradians 

As  above. 

Energies  Available 
(Single  shot  only) 

Laser  Measured 

Setting  Values 

0.15  0.14 

0.25  0.34 

1.0  0.43 

2.5  0.74 

4.0  1.04 

5.0  1.96 

7.5  2.44 

10.0  2.09 

All  values  listed  in 
millijoules. 

Determined  experimentally 
for  532  nanometers  using  a 

1 .06  micron  laser  safety 
filter. 
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4  mm 


Off-axis  artifact  approximately  1  mm  in 
diameter 

Figure  6.  Nd:YAG  Laser  Beam  profile. 


2 .  Dall-Kirkham  Telescope 

The  Dall-Kirkham  telescope  is  a  form  of  Cassegrain  telescope  but  the  primary 
mirror  is  aspheric  and  the  secondary  is  spherical.  A  summary  of  the  telescope’s 
specifications  is  in  Table  2. 

The  telescope  was  used  initially  for  atmospheric  studies  in  the  1970s  and  it  was 
refurbished  for  use  in  this  lidar  applicadon.  The  original  mirror  smoothness  was  quoted  as 
a  tenth  of  a  wavelength  at  632.8  nanometers  but  an  inspection  of  the  mirrors  found  a 
number  of  pits  in  the  coating.  It  is  probable  that  the  smoothness  is  no  longer  to  this 
precision  due  to  salt  corrosion;  but  for  this  lidar  applicadon  the  existing  mirror  quality  will 
be  more  than  adequate  since  the  telescope  is  not  being  used  for  imaging.  Also  any 
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additional  mirror  losses  due  to  the  corrosion  can  he  made  up  easily  by  increasing  the  laser 
output  energy.  The  telescope  has  a  4  inch  diagonal  to  redirect  the  focal  plane  90  degrees  to 
the  telescope  axis  and  it  will  he  also  used  for  the  lidar  application. 


TABLE  2.  DALL-KIRKHAM  TELESCOPE  SPECIFICATIONS  [Ref.  28) 


Manufacturer 

Tinsley  Laboratories  Incorporated 

3900  Lakeside  Avenue 

Richmond,  California  94806 

In  approximately  early  1970s 

Primary  Mirror 

Diameter  -  17.75  inches  (aspheric) 

Radius  of  Curvature  -  1 26  inches 

Focal  length  -  63.0  inches 

Optical  coating  -  Al/SiO 

Substrate  -  Cervit 

F  number  -  3.55 

Secondary  Mirror 

Diameter  -  5.75  inches  (spherical) 

Radius  of  Curvature  -  44.80  inches 

Optical  coating  -  Al/SiO 

Substrate  -  Cervit 

Overall  System 

Effective  focal  length  -  252  inches 

Mirror  separation  -  46.25  inches 
(vertex  to  vertex) 

Effective  F  number  -  15.0 

Focal  plane  is  21.0  inches  behind  vertex  of 
primary  mirror. 

The  telescope  was  mounted  on  a  Pelco  PT2000L  Pan  and  Tilt  Unit  to  have  the 
ability  to  align  the  lidar  on  targets  of  interest.  The  unit  has  a  ±  90°  adjustable  tilt  and  355° 
adjustable  pan  movement  [Ref.  29).  A  camera  and  surplus  right-angle  telescope  were 
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available  for  mounting  on  to  the  telescope  to  provide  an  indication  of  the  telescope’s  field  of 
view  as  well  as  a  warning  of  obstructions  but  time  did  net  permit  their  installation. 

3 .  Laser  Optics 

The  lidar  transmitter  laser  optics  layout  is  sho' vn  at  Figure  7.  The  1 .06  micron 
wavelength  beam  was  split  from  the  532  nanometer  main  laser  beam  and  directed  into  a 
beam  stop  using  a  1.06  micron  reflecting  beam  splitter.  The  provision  was  made  for  future 
use  of  the  1.06  micron  wavelength,  if  required,  by  the  addition  of  infrared  optics.  The 
choice  of  the  beam  expanding  optics  was  very  important  in  order  to  properly  fill  the 
telescope’s  secondary  mirror  without  loss  of  energy.  In  the  calculations  for  the  beam 
expander  optics,  it  was  decided  to  transmit  as  much  laser  energy  as  possible  into  the 
atmosphere.  Therefore,  the  99%  beam  radius  point  was  taken  as  pi  times  beam  waist 
divided  by  two  (  or  3.141  millimeters).  The  radius  of  the  secondary  mirror  was  chosen  as 
the  99%  point  and  the  distance  to  focal  point  was  used  in  the  calculations.  The  matrix 
representation  for  the  solution  using  only  one  lens  is  shown  below: 
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w’  is  the  beam  radius  on  the  secondary  mirror  at  the  1/e  point  (46.49  millimeters), 

0’  is  the  beam  half  angle  beam  divergence  needed  to  fill  the  secondary  at  the  1/e  point 
and  knowing  the  focal  plane  location  (27.229  milliradians), 

T  is  the  distance  from  the  secondary  mirror  to  lens, 
f  is  the  focal  length  of  lens. 

w0  is  the  measured  laser  beam  waist  (radius  2.0  millimeters),  and 
0O  is  the  measured  laser  half  angle  beam  divergence  (0.165  milliradians). 
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To  Telescope 

^  Opening  through  Box  Wall 


Figure  7.  Laser  Box  Optics. 


Trigger 


The  lenses  used  in  the  beam  expander  were  restricted  to  those  readily  available 
within  the  NPS  Department  of  Physics.  A  -74  millimeter  focal  length  meniscus  lens  was 
dete-mined  to  be  the  best  design  to  expand  the  laser  beam  from  Equation  (3-1).  The 
distance  for  the  lens  from  the  secondary  mirror  was  determined  to  be  163.4  centimeters 
using  Equation  (3-1).  Other  designs  using  multiple  lenses  were  considered  but  disregarded 
as  showing  no  advantage  over  the  simpler  single  lens  approach. 
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The  Nd:YAG  laser  had  the  capability  to  incorporate  a  helium  neon  laser  as  a 
pointing  device;  this  capability  consisted  of  a  electrical  power  connection  and  the  front  and 
rear  Nd:YAG  mirrors  were  transmissive  at  632.8  nanometers.  The  helium  neon  laser  was 
not  purchased  with  the  Nd:YAG  laser,  so  a  visible  670  nanometer  red  laser  diode  was 
substituted.  The  laser  diode  was  used  as  a  diagnostic  tool  for  the  alignment  of  die  NdiYAG 
laser  beam  in  the  laser  box  without  constantly  firing  the  Nd:YAG  laser.  A  9  volt  battery 
was  used  to  power  the  laser  diode  in  order  to  have  an  independent  power  supply.  The  9 
volt  power  supply,  without  the  battery,  was  left  inside  the  sealed  laser  box  for  future  beam 
alignment.  There  was  a  concern  about  damage  to  the  laser  diode  due  to  the  1.06  micron 
laser  beam  energy  leakage  from  the  Nd:YAG  rear  mirror,  which  was  measured  at  a 
constant  0.25  millijoules  for  all  single  pulse  shot  settings.  A  helium  neon  laser  mirror, 
transmissive  at  1.06  microns,  was  placed  directly  behind  the  rear  mirror  of  the  Nd:YAG 
laser,  as  shown  in  Figure  7,  to  prevent  the  1.06  micron  laser  beam  energy  leakage  from 
damaging  the  laser  diode. 

A  surplus  PIN  photodiode  was  added  into  the  design  in  order  to  trigger  the 
digital  oscilloscope  for  lidar  backscatter  measurements.  The  PIN  photodiode  was  placed 
inside  the  laser  box  and  it  was  aimed  at  one  of  the  mirrors.  There  was  enough  energy 
scattered  from  the  mirrors  for  the  PIN  diode  to  produce  a  signal  which  would  trigger  the 
oscilloscope.  The  Nd:YAG  Hash  lamp  was  observed  at  the  PIN  photodiode  output  and  a 
highpass  filter  was  built  using  a  floating  ground  to  isolate  the  Nd:YAG  laser  pulses.  The 
floating  ground  was  required  because  the  very  large  current  surge  from  the  flash  lamp 
caused  an  induced  surge  in  the  common  ground.  The  PIN  photodiode  and  associated 
highpass  filter  were  then  used  to  trigger  the  digital  oscilloscope  to  start  the  lidar 
measurements. 
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After  initial  testing  of  the  laser  box  mounted  with  the  telescope,  an  overfilling 
of  the  secondary  mirror  was  observed  which  in  the  far-lleld  appeared  as  an  expanding  ring 
encircling  the  collimated  lidar  beam.  This  overfilling  was  attributed  to  a  slight 
misalignment  of  the  beam  expanding  lens  (which  was  later  corrected)  and  the  non-Gaussian 
laser  beam  shape  causing  the  laser  beam  edges  to  be  of  a  higher  energy  density  than 
expected.  It  was  decided  to  employ  a  baffle,  or  iris,  to  restrict  the  exit  aperture  out  of  the 
laser  box  to  that  needed  to  fill  the  secondary  mirror. 

The  half-angle  beam  divergence  for  the  lidar  transmitter  was  determined  to  be 
0.794  milliradians  by  measuring  the  growth  of  the  laser  beam’s  diameter  over  a  known 
distance.  The  beam  divergence  was  decided  to  be  adequate  for  lidar  use  and  safety 
considerations,  since  the  laser  beam  would  never  focus  in  the  far-field. 

C.  LASER  SAFETY 

Laser  safety  was  an  important  consideration  in  the  design  of  the  lidar.  In  the 
designed  configuration,  the  laser  beam  must  be  greatly  expanded  using  a  telescope  in  order 
to  achieve  a  maximum  allowed  energy  density,  the  Maximum  Permissible  Exposure 
(MPE),  of  5  x  10'7  joules  per  centimeter2  [Ref.  24].  The  energy  density  of  the  transmitted 
laser  beam  was  measured  using  a  Laser  Precision  Corporation  energy  meter  which  had  a 
200  microjoule  range  and  a  five  square  centimeter  detection  area.  Since  the  secondary 
mirror  blocked  the  transmission  of  the  peak  laser  energy  along  the  optical  axis,  the  peak 
transmitted  energy  density  was  found  along  the  edge  of  the  secondary  mirror  as  expected. 
The  measured  peak  energy  densities  in  the  expanded  laser  beam  for  different  laser  setting  is 
shown  in  Table  3. 

The  range  along  the  unobstructed  optical  axis  of  the  laser  beam,  wb'ch  it  can  be 
safely  viewed  by  the  human  eye,  is  referred  to  as  the  Nominal  Ocular  Hazard  Distance 
(NOHD)  [Ref.  24].  The  NOHD  defines  the  range  at  which  the  wavelength  specific  MPE  is 
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not  exceeded.  This  distance  is  one  of  the  principle  measures  of  laser  eye  safety  and  it  will 
be  determined  for  each  laser  setting.  The  NOHD  can  be  calculated  from  the  average  radiant 
exposure.  Equation  B9  in  reference  24,  which  is  defined  as  the  energy  of  the  laser  beam 
divided  by  its  total  area: 


H  =  0  exp  [-^r] 
[(a  +  r  0) 


(3-3) 


where 


H  is  the  radiant  exposure  measured  in  joules  per  centimeter^, 

Q  is  the  radiant  energy  output  of  a  pulsed  laser  measured  in  joules 
(defined  at  the  1/e  points  of  maximum  irradiance), 

(i  is  the  atmospheric  extinction  coefficient, 
r  is  the  range  in  centimeters, 

a  is  the  laser  beam  diameter  measured  to  1/e  points  in  centimeters,  and 
<J)  is  the  laser  beam  divergence 

(twice  the  half-angle  beam  divergence,  1 .588  milliradians). 


In  the  case  of  this  application,  the  lidar  has  a  secondary  mirror  which  blocks  the 
center  transmission  of  laser  beam,  as  shown  in  Figure  5.  The  total  radiant  energy  output, 
Q,  which  actually  exits  the  lidar  transmitter,  can  be  determined  from  the  fraction  of  the  total 
power  which  is  blocked  by  the  secondary  mirror  and  noting  that  Q  is  equal  to  total  power 
divided  by  laser  pulse  duration.  The  fraction  of  the  total  power  or  energy  of  a  Gaussian 
beam,  which  is  blocked  by  the  secondary  mirror,  is  obtained  using  the  following  relation: 


fd  =  l  -  exp  =  0.4045  (3-4) 

F0  \  Dl2  / 
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where 


Pd  is  the  power  blocked  by  the  secondary  mirror, 

P0  is  the  total  power  of  the  laser  beam, 

D0  is  the  diameter  of  the  secondary  mirror  ( 14.61  centimeters),  and 
Dl  is  the  laser  beam  diameter  measured  from  the  1/e  points 

(28.7  centimeters  -  twice  the  diameter  of  telescope  divided  by  n , 
since  the  telescope  was  designed  to  transmit  99%  of  the  laser  energy). 


The  denominator  of  Equation  (3-3)  must  also  be  modified  to  account  for  the 
reduction  of  laser  beam  area  due  to  the  blocking  of  the  secondary  mirror.  The  modified 
laser  beam  area  can  be  simply  written,  using  previously  defined  terms,  as: 


K 


(a  +  r  <p) 
2 


iz 


(Dl  +  r  <J>) 
2 


(3-5) 


An  expression  for  the  lidar  radiant  exposure,  //,  can  be  written  by  substituting 
Equation  (3-5),  and  the  value  from  Equation  (3-4)  into  Equation  (3-3).  Q  can  be  written  in 
terms  of  the  total  laser  beam  energy  by  noting  that  86%  of  the  laser’s  energy  is  found 
between  the  1/e  points  of  maximum  irradiance;  Q  =  0.86  times  the  total  laser  beam  energy. 
The  resulting  new  expression  for  H  is: 


( 1  -  0.4045)  0.86  Qt  exp  [-[ir] 


2 

n 

(Dl  +  r<f>) 

-  7t 

D* 

2 

L  2  J 

2  J 

(3-6) 


where  Qr  is  the  total  energy  of  the  laser  beam. 
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TABLE  3.  LASER  SAFETY 


Laser  Setting 

Total  Laser  Energy 
Measured  at  Laser 
Box  Exit  -  5  shot 
average  (see  text). 

Or 

(milli  joules) 

Peak  Energy 
Density  Measured 
using  5  cm2  Energy 
Probe  across 
Telescope. 

(microjoules  per 
cm2) 

NOHD 

(meters) 

Worst  Case  -  no 
atmospheric 
attenuation. 

(see  text) 

10 

2.414 

2.0 

184.1 

7.5 

2.06 

2.96 

158.1 

5.0 

1.96 

1.58 

150.4 

4.0 

1.123 

1.34 

77 

2.5 

1.7 

1.125 

129.7 

1.8 

1.885 

0.6725 

144.5 

1.2 

1.61 

121.9 

0.9 

0.917 

wgggm 

55.5 

A  worst  case  estimate  for  the  NOHD  can  be  obtained  from  Equation  (3-6)  if  the 
atmospheric  extinction  coefficient  term  is  neglected,  exp  [-pr],  and  H  is  set  to  the  MPE. 
The  resulting  first  order  differential  equation,  in  terms  of  the  range,  can  be  solved  for  each 
total  laser  energy  output  as  shown  in  Table  3: 
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(3-7) 


r,  )  2Dl  r  |  (d£-  Dp)  (1  -  0.4045)  0.86  Qr 
4>  (j,2  JMPE 

Upon  closer  inspection  of  Table  3.  there  were  observed  a  number  of  inconsistencies 
in  the  measured  energy  densities  for  laser  settings  0.9  and  1 .2  and  the  calculated  NOHD. 
The  values  for  the  NOHD  should  have  been  very  small  based  on  the  MPE,  but  they  were 
not.  A  more  realistic  indication  of  laser  safety  would  the  peak  energy  density  measured 
across  the  telescope  as  shown  in  the  third  column  of  Table  3. 

It  was  noted  that  the  measured  laser  box  energy  values,  shown  in  the  second  column 
of  Table  3,  were  different  from  those  measured  in  Table  1.  This  could  be  attributed  to 
possible  misalignment  of  the  neuual  density  filters  which  are  mounted  in  the  Nd:YAG  laser 
and  provide  the  selectable  energy  settings.  This  problem  could  not  be  repaired  by  the 
author. 

D.  LIDAR  RECEIVER 

The  designed  lidar  receiver  layout  is  at  Figure  8.  The  receiver  was  designed  to  have 
a  larger  field  of  view  than  the  transmitter  to  ensure  that  as  much  of  the  backscattered  photon 
flux  as  possible  would  be  collected.  The  receiver’s  collecting  lens  had  a  focal  length  250 
millimeters  and  had  a  diameter  of  145  millimeters. 

There  were  two  techniques  employed  for  reducing  the  amount  of  background  sky 
radiance  collected  by  the  receiver,  a  pin-hole  and  a  narrowband  filter.  An  one-eighth  inch 
diameter  pin-hole  was  placed  in  the  focal  plane  of  the  receiver’s  collecting  lens.  The 
resulting  reduced  field  of  view  was  calculated  to  be  0.127  x  10-3  steradians  or  a  half  angle 
subtense  of  6.35  milliradians.  A  smaller  pin-hole  could  be  used  to  improve  the  far-field 
lidar  signal  by  reducing  the  background  in  the  field  of  view  and  thereby  improving  the 
SNR.  A  ten-nanometer  narrowband  optical  filter  centered  at  532  nanometers  was  also  used 
in  the  lidar  receiver,  which  reduced  the  overall  background  radiance  to  a  very  narrow  band 


38 


around  the  laser  wavelength.  A  collimating  lens  was  positioned  in  front  of  the  narrowband 
filter  to  ensure  that  all  of  the  collected  photons  passing  through  the  pin-hole  would  strike 
the  photocathode  and  not  miss  it  due  to  divergence. 

The  photomultiplier  tube  used  was  a  Hamamatsu  R1913.  It  was  chosen  for  its  fast 
rise  response  time  and  its  sensitivity  in  the  visible  region.  The  specifications  of  the  R 19 13 
photomultiplier  tube  are  in  Table  4. 


Figure  8.  Lidar  Receiver. 

The  output  of  the  photomultiplier  tube  was  to  be  amplified  using  a  Model  VV100B 
wideband  pulse  amplifier  made  by  the  LeCroy  Corporation  in  the  initial  design  of  the  lidar 
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receiver.  The  amplifier  has  a  gain  of  10,  less  than  2  nanosecond  rise  time  and  >  200 
megaHertz  bandwidth  [Ref.  30).  The  LeCroy  Corporation  also  makes  a  mounted  version 
of  the  Model  VV100B  (which  has  the  designation  Model  VV100BTB)  and  it  had  the  proper 
ground  planes  and  components  for  low  noise  operation  which  made  it  ideal  as  a  lidar  signal 
amplifier.  The  amplifier  was  not  required  in  the  final  design  because  the  digital 
oscilloscope  used  to  acquire  the  data  was  very  sensitive;  also  there  was  a  lot  of  noise  in  the 
lidar  system.  The  source  of  the  electrical  noise  was  not  found  and  it  was  decided  that  the 
amplifier  would  only  compound  the  problem  by  amplifying  the  noise.  The  amplifier  is 
planned  for  use  in  a  modified  design  of  this  lidar  where  photon  counting  lidar  techniques 
are  to  be  used  after  the  noise  source  is  eliminated. 

TABLE  4.  R1913  PHOTOMULTIPLIER  SPECIFICATIONS  [Refs.  31  and  32] 


Cathode  Sensitivity 

65  miiliAmperes  per  Watt  (typical) 

55  miiliAmperes  per  Watt  (532  nanometers) 

Current  Gain 

3.0  x  103 

Response  Time 

1.0  nanoseconds  (at  2000  VDC) 

Anode  Dark  Current 

Typical  0.05  nanoAmperes 

Maximum  0.5  nanoAmperes 

Spectral  Response 

Range  185-900  nanometers 

Peak  Wavelength  420  nanometers 

Quantum  Efficiency 

Average  in  the  Visible  Region  13% 

At  532  nanometers  15  % 
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E.  LIDAR  ALIGNMENT 


1 .  Telescope 

The  telescope  was  collimated  for  the  lidar  application,  employing  an  amateur 
astronomer  technique  [Ref.  33]  without  the  4  inch  diagonal  mirror,  using  a  sufficiently 
bright  star  in  the  Monterey  Bay  .sky  to  obtain  an  observable  focused  point  of  light  at  the 
focal  plane.  The  focal  plane  was  then  adjusted  by  moving  the  secondary  mirror  so  that  the 
focal  plane  was  21.0  inches  behind  the  vertex  of  the  primary  mirror  as  given  in  the 
telescope’s  specifications.  The  mechanical  alignment  of  all  the  mirrors  was  accomplished 
by  observing  the  double  reflection  from  the  secondary  and  primary  mirrors  of  the  author’s 
eye  centered  in  the  diagonal  mirror  . 

The  confirmation  of  the  collimation  was  done  by  using  a  second  smaller 
telescope,  which  was  focused  at  infinity.  A  light  source  was  placed  in  the  large  telescope’s 
focal  plane  and  then  the  smaller  telescope  was  placed  in  the  large  telescope’s  field  of  view. 
The  theory  for  this  method  of  collimation  was  that  if  the  light  source  was  at  the  focal  plane 
of  the  large  telescope,  then  the  image  of  the  light  formed  by  the  smaller  telescope  would  be 
in  focus,  and  it  was. 

2.  Lidar  Transmitter  and  Receiver  Alignment 

The  lidar  transmitter  was  initially  aligned  using  the  red  laser  diode  shown  in 
Figure  7.  In  theory,  if  the  red  laser  diode  was  aligned  with  the  Nd:YAG  laser,  the  mirrors 
could  be  adjusted  so  that  the  unexpanded  red  laser  diode  reflection  off  the  secondary  mirror 
followed  its  path  back  into  the  laser  box.  Then,  the  Nd:YAG  laser  axis  was  aligned  to  the 
telescope’s  secondary  mirror  axis.  In  practice,  the  alignment  of  the  red  laser  diode  to  the 
Nd:YAG  laser  axis  was  never  really  achieved  because  of  problems  with  precisely  mounting 
the  laser  diode.  The  laser  diode  was  able  to  provide  a  reasonably  accurate  pointing  for  the 
Nd:YAG  laser  and  this  was  used  to  roughly  align  the  Nd:YAG  optical  axis  with  the 
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telescope’s  optical  axis.  A  more  precise  alignment  was  conducted  by  measuring  the  energy 
density  distribution  around  the  secondary  mirror  obstruction.  If  the  lidar  transmitter  was 
aligned,  then  the  energy  density  would  follow  the  energy  pattern  in  Figure  6.  A  baffle  was 
used  to  restrict  the  size  of  the  expanded  Nd.YAG  laser  beam  as  mentioned  in  Section  B. 
The  telescope’s  secondary  mirror  was  slightly  adjusted  to  compensate  for  any  inaccuracies 
in  the  laser  box  alignment.  The  telescope’s  primary  mirror  was  never  adjusted.  The  beam 
shape  and  energy  distribution  of  the  collimated  laser  beam  was  observed  by  using  a  locally 
made  cardboard  screen  painted  florescent  yellow.  The  yellow  florescence  after  irradiation 
by  the  laser  beam  was  easily  observed  through  the  laser  safety  goggles  and  this  was  a  very 
inexpensive  method  to  view  the  beam. 

The  lidar  receiver  was  aligned  after  the  lidar  transmitter  alignment.  It  was 
accomplished  by  adjusting  the  position  of  the  receiver  pin-hole  to  obtain  the  largest  returned 
signal  from  the  painted  cardboard  screen  in  the  laboratory.  The  procedure  was  again 
repeated  once  the  lidar  was  installed  on  the  roof  of  Spanagel  Hall  using  the  same  painted 
cardboard  screen  but  at  a  longer  distance. 

F .  DATA  ACQUISITION 

The  measurement  and  digitizing  of  the  backscattered  lidar  signal  was  accomplished 
using  a  Tektronix  DSA  602A  Digitizing  Signal  Analyzer  which  had  a  two  gigasample  per 
second  sampling  rate  [Ref.  34],  A  signal  bandwidth  of  1  gigaHertz  was  obtained  by  using 
an  1 1 A72  two-channel  plug-in  module.  The  Tektronix  DSA  602A  had  a  MS-DOS  3.5  inch 
disc  drive  built  into  it  which  was  a  useful  feature.  This  allowed  data  storage  directly  to 
disc,  otherwise  the  data  would  have  to  be  transferred  to  a  computer  and  then  to  disc  which 
would  be  a  slower  process.  The  range  resolution  for  DSA  602A  was  calculated  to  be  0.15 
meters  but  this  is  the  two-way  range  resolution  and  needs  to  be  divided  by  two  to  reflect  the 
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range  at  which  the  backscattered  signal  originates.  The  actual  timebase  range  resolution  is 
dominated  by  the  laser  pulse  width  (t)  and  is  estimated  as  0.6  meters  (cx/2). 

The  Tektronix  DSA  602A  became  available  at  the  last  minute,  but,  initially  it  was 
planned  to  use  the  Hewlett  Packard  HP54111D  Digitizing  Oscilloscope  which  had  a 
sampling  rate  of  1  gigasample  per  second.  In  order  to  store  the  data  for  off-line  analysis 
using  the  Hewlett  Packard  HP541 1  ID,  a  computer  interface  was  required  as  described  in 
the  previous  paragraph.  The  Tektronix  DSA  602A  was  chosen  over  the  Hewlett  Packard 
HP54111D  mainly  because  of  the  ease  of  data  storage.  The  theoretical  lidar  inversion 
analysis  conducted,  in  Appendix  B,  used  the  Hewlett  Packard  HP54111D  sampling  rate 
which  provided  a  one-way  timebase  range  resolution  for  the  lidar  inversion  of  0.15  meters. 
It  was  decided  not  to  change  the  analysis  to  reflect  the  new  range  resolution  since  the  small 
change  would  not  have  produced  any  better  results. 

The  lidar  backscattered  signal  could  have  been  observed  more  easily  if  a  logarithmic 
amplifier  was  used  at  the  oscilloscope  input.  The  logarithmic  amplifier  would  expand  the 
lidar  backscattered  display  by  amplifying  the  far-field  return  while  leaving  the  near-field 
return  in  its  original  form.  This  would  allow  both  the  near-  and  far-field  returns  to  be 
observed.  The  data  could  be  reverted  to  its  original  form  by  using  the  manufacturer’s 
logarithmic  amplifier  specifications.  This  would  be  a  possible  future  design  modification. 

The  next  chapter  will  present  various  actual  measured  lidar  signals  from  hard  targets 
and  clouds. 
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IV.  LIDAR  DEMONSTRATION  AND  RESULTS 

A.  DEMONSTRATION  OVERVIEW 

On  9  June  1993  at  1945  hours,  the  lidar  designed  in  Chapter  HI  was  tested  to  ensure 
that  it  could  provide  range  data  to  clouds  which  was  the  primary  aim  of  this  thesis.  The 
weather  conditions  were  clear  with  light  clouds  and  fog  over  Monterey  Bay.  There  was  a 
fog  bank  forming  off  the  coast,  and  later  in  the  evening  it  was  ranged  to.  The  lidar  was 
operated  under  strict  control  of  the  NPS  Physics  Department  staff  and  the  author  was 
present  to  take  data.  The  data  was  recorded  using  the  Tektronix  DSA  602A  Digitizing 
Signal  Analyzer  as  described  in  Chapter  III.  Both  the  lidar  receiver  and  transmitter  trigger 
signals  were  stored  by  the  oscilloscope  and  the  trigger  signal  was  used  in  the  later  analysis 
to  confirm  the  start  time  for  the  lidar  return.  The  aim  of  the  demonstration  was  to  obtain  a 
range  to  hard  targets  initially  and  then  to  clouds. 

B .  RANGE  RESULTS 

The  lidar  was  initially  tested  against  hard  targets  to  determine  if  a  reflected  signal 
could  be  observed  by  the  lidar  receiver.  A  tree  on  the  school  grounds  was  used  as  the  first 
hard  target,  estimated  at  about  300  meters  from  the  lidar.  The  lidar  transmitter  range 
domain  signal  is  displayed  in  Figure  9  for  the  tree  target.  The  tree  return  is  clearly  visible 
in  Figure  9  and  this  confirmed  the  lidar’s  ranging  ability  for  short  ranges.  The  digital 
oscilloscope  was  not  used  in  its  highest  resolution  because  this  was  a  proof  of  concept 
demonstration  only  vice  actual  atmospheric  studies.  Ranges  were  also  obtained  to  the 
Jacks  Peak  Ridge  which  runs  southwest  of  the  city  of  Monterey.  The  range  measured  from 
a  city  of  Monterey  map  indicate  the  range  from  the  school  to  the  ridge  to  be  in  the  order  of 
five  kilometers.  The  ranges  observed  were  approximately  six  kilometers  which  are 
considered  more  accurate. 
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Figure  9.  Lidar  return  from  tree  on  campus  grounds  at  approximately  2015  hours. 

The  next  demonstration  of  the  lidar  was  the  ranging  to  clouds.  During  the  course  of 
the  evening,  a  number  of  light  cumulus  clouds  began  to  blow  in  from  the  ocean  and  they 
were  successfully  ranged.  Figure  10  is  an  example  of  one  such  cloud  ranging.  This  figure 
also  has  some  interesting  features  such  as  the  characteristic  near-field  exponential  decay  of 
the  lidar  signal  as  it  propagates  away  from  the  lidar  receiver  and  a  strong  cloud  lidar  signal. 
The  cloud  was  visually  estimated  to  be  one  kilometer  and  the  actual  range  determined  from 
Figure  10  was  1400  meters.  The  lidar  was  able  to  successfully  range  to  clouds  with  good 
agreement  from  visual  estimation.  A  number  of  other  clouds  were  ranged  during  the 
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Lidar  signal  (volts) 


course  of  the  evening  and  a  variety  of  lidar  returns  were  observed.  There  were  variations 
in  the  magnitudes  of  the  lidar  cloud  returns  which  were  attributed  to  the  changing  nature  of 
the  atmosphere  due  to  a  strong  breeze  blowing  in  from  the  ocean. 


Range  (meters) 

Figure  10.  Lidar  return  from  a  cloud  at  approximately  2040  hours. 

The  last  lidar  measurements,  taken  at  approximately  2110  hours,  were  of  the  fog 
bank  off  the  coast  which  was  observed  forming  earlier.  The  fog  bank  lidar  return  is  shown 
in  Figure  11. 
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Figure  LI.  Lidar  return  from  fog  bank  at  approximately  2110  hours. 

It  was  noted  that,  as  the  evening  progressed  and  the  sun  set,  the  background  noise 
decreased  as  expected,  which  can  be  observed  by  comparing  of  the  noise  levels  in  Figures 
9  and  11.  The  noise  floor  in  Figure  10  was  further  reduced  by  averaging  two  consecutive 
lidar  returns  from  the  same  cloud.  In  order  to  conduct  daytime  lidar  analysis,  the  noise 
levels  will  have  to  be  reduced  by  employing  a  smaller  pin-hole  as  discussed  in  Chapter  III 
as  well  as  sealing  the  lidar  receiver  enclosure  better  to  prevent  stray  light  from  entering  the 
lidar  receiver  box.  An  averaging  technique,  such  as  photon  counting,  could  be  employed 
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to  reduce  the  daytime  noise  levels  but  typically  the  pulse  repetition  rates  are  in  the  kiloHertz 
range  to  produce  good  results. 

A  further  discussion  of  the  theoretical  day  and  nighttime  lidar  ranges  is  now 
necessary  for  comparison  purposes.  The  daytime  maximum  theoretical  range,  as  derived  in 
Equauon  (2-27),  can  be  determined  by  substitudng  in  the  design  and  measured  values  from 
Chapter  HI: 


(4-1) 


where 


Tnuw  is  the  maximum  range  for  the  lidar  system  (meters), 

PQ  is  the  lidar  system  transmitted  power 

(0.5  x  106  watts,  2  millijoules/4  nanoseconds), 

T  is  the  laser  pulse  duradon  (4  nanoseconds), 

B(rmax)  is  the  atmospheric  volume  backscattering  coefficient 

(a  range  of  values  are  given  in  reference  25  between  1  x  1H7  and  1  x  10-4 
meters1  •  steradians1.  A  value  of  6  x  10-6  meters1  •  steradians1 
will  be  assumed), 

c  is  the  speed  of  light  (3  x  10*  meters/second), 

Aq  is  the  receiver’s  effecuve  aperture  (0.0165  meter2), 

X  is  the  laser  wavelength  (532  x  10-9  meters), 
tj(X)  is  the  photomuldplier  tube  quantum  efficiency  (15%), 
a(r)  is  the  atmospheric  volume  exuncdon  coefficient  (assumed 
to  be  constant  along  the  path  at  2  x  10-3  meters-1  [Ref.  35]), 
h  is  Planck’s  constant  (6.6266176  x  IO-34  watts  •  seconds2), 

AX  is  the  narrowband  filter  width  (0.010  microns), 

Sb(X)  is  the  background  spectral  radiance  (10  watts  meter-2  •  steradians1  •  micron-1 
[Ref.  22]), 

To(X)  is  the  receiver  opdcs  transmittance  (assumed  to  be  0.8), 

Qo  is  the  receiver  solid  angle  (0. 127  x  I03  steradians),  and 
Af  is  the  electrical  system  bandwidth  (assumed  to  be  1  gigaHertz). 
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The  resulting  expression  for  the  maximum  nighttime  range  is: 

rmax  =  1.0137  x  HI3  exp  [-2  x  10  3  rmax]  (meters).  (4-3) 

The  nighttime  maximum  theoretical  range,  as  derived  in  Equation  (2-32),  can  be 
determined  by  substituting  in  the  design  and  measured  values  from  Chapter  HI: 

(4-4) 

The  resulting  expression  for  the  maximum  nighttime  range  is: 

rmax  =  6.305  x  lO3  exp  [-2  x  10-3  rmaxJ  (meters).  (4-5) 

Equations  (4-3)  and  (4-5)  must  be  solved  iteratively  to  obtain  a  value  for  rmax.  The 
daytime  and  nighttime  maximum  lidar  ranges  are  0.9  and  3.3  kilometers  respectively.  The 
nighttime  maximum  range  seems  to  be  low  when  comparing  it  to  the  lidar  return  in  Figure 
1 1.  This  should  not  be  possible  as  the  lidar  appears  to  outperform  its  theoretical  limits. 
Therefore,  the  overall  system  bandwidth  may  be  less  than  1  gigaHertz  or  the  magnitudes  of 
P(r)  and  o(r)  may  have  been  incorrect  for  the  weather  conditions.  The  values  chosen  for 
P(r)  and  o(r)  reflect  the  average  values  for  clear  air  conditions  but  the  conditions  may  have 
been  better  than  estimated  in  the  near- field. 


C.  LIDAR  INVERSION 

A  lidar  inversion  was  carried  out  on  the  lidar  return  shown  in  Figure  10  as  it  was  one 
of  the  best  lidar  returns.  The  algorithms  developed  in  Chapter  II  and  Appendix  B  were 
used  to  conduct  the  lidar  inversion.  It  should  be  noted  that  the  range  resolution  for  the  lidar 
return  was  1.2  meters  which  was  well  above  the  high  resolution  theoretical  lidar  inversions 
carried  out  in  Appendix  B.  Therefore,  it  will  not  be  the  aim  to  produce  an  accurate  lidar 
inversion  but  to  demonstrate  the  capability. 

1 .  Determination  of  r„ 

Before  the  lidar  inversion  can  progress  a  value  for  rD,  the  receiver/transmitter 
overlap  range,  needs  to  be  determined.  This  was  done  theoretically  by  measuring  the 
displacement  of  the  receiver  and  transmitter  axis,  the  pin-hole  radius  in  the  receiver  which 
will  provide  the  field  of  view  and  the  transmitter  half-angle  beam  divt  rgence.  From 
Chapter  III,  the  pin-hole  diameter  was  set  at  1/8”  and  the  transmitter  half  angle  beam 
divergence  was  measured  at  0.794  milliradians.  The  displacement  of  the  two  axes  was 
35.56  centimeters.  The  graphical  solution  for  r0  showing  both  the  receiver’s  fields  of  view 
and  transmitter’s  beam  expansion,  in  two  dimensions,  is  in  Figure  12. 

Therefore,  the  minimum  theoretical  range  at  which  the  receiver’s  field  of  view 
overlaps  the  transmitter’s  beam  is  104.6  meters.  This  value  will  be  used  in  the  lidar 
inversion  process. 
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Beam  Size  (meters) 


2.  Lidar  inversion 

The  inversion  of  the  lidar  signal  is  not  as  straightforward  as  presented  in 
Chapter  II  and  Appendix  B.  The  determination  of  the  lidar  constant,  C ,  from  Equation  (2- 
10)  can  be  difficult  without  prior  knowledge  for  the  extinction  coefficient  profile.  If  the 
extinction  coefficient  is  known,  then  the  Lidar  equation  can  be  worked  backwards  to 
determine  the  lidar  constant  which  can  be  used  for  future  lidar  inversions.  The  receiver 
output  signal  is  a  function  of  the  incident  power,  as  presented  in  Equation  (2-25),  but  in 
reality  there  is  still  an  overall  system  efficiency  which  must  be  considered.  Therefore,  a 
practical  solution  to  the  problem  of  the  lidar  constant  and  receiver  efficiency  would  be  to 
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group  these  constants  together  and  determine  them  as  a  whole  by  calibrating  the  lidar 
against  a  known  extinction  coefficient  profile.  The  only  problem  with  this  approach  would 
be  acquiring  the  necessary  equipment,  such  as  a  visioceilometer  as  described  in  reference  9. 

For  the  purpose  of  this  lidar  inversion,  the  lidar  constant  will  be  set  to  an 
arbitrary  value.  The  lidar  inversion  will  be  inaccurate  at  the  boundary  but  the  large  range 
resolution  will  make  an  accurate  lidar  inversion  impractical  anyway.  The  inverted  lidar 
signal  is  shown  at  Figure  13. 


range-ro  (km) 

Figure  13.  Lidar  inversion  of  Figure  1 1 
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The  lidar  signal,  Irom  Figure  1 1,  was  first  smoothed  by  piecewise  curve  fitting 
to  remove  the  noise  from  the  lidar  signal.  The  cloud  return,  shown  in  the  center  of  Figure 
1 1 ,  was  left  to  get  a  comparison  of  smoothed  and  raw  lidar  data  inversion.  The  piecewise 
smoothing  process  was  time  consuming  and  investigations  into  a  more  practical  method  are 
needed.  The  extinction  coefficient  profile  is  shown  in  Figure  13  and  it  became  very  large  at 
the  boundary.  This  is  a  concern  because  Figure  1 1  does  not  show  any  returns  past  the 
initial  cloud  return  at  1300  meters.  Therefore,  this  is  the  boundary  inaccuracy  mentioned 
earlier  due  to  the  arbitrary  selection  of  the  lidar  constant  and  the  large  range  resolution. 

The  cloud  irradiated  by  the  lidar  in  Figure  1 1  was  a  light  cumulus  cloud.  The 
lidar  inversion  of  the  cloud  region  is  reasonable  since  the  expected  extinction  coefficient  in 
cumulus  clouds  is  in  the  order  of  10  kilometer-1  [Ref.  36].  The  lidar  inversion  peak 
measured  extinction  coefficient  was  5  kilometer-1  which  was  close  to  the  expected  value. 
There  were  errors  in  the  lidar  inversion  as  mentioned  earlier  but  the  profile  of  the  cloud  was 
preserved  to  demonstrate  the  robustness  of  the  process. 

The  designed  lidar  met  its  objective  to  provide  accurate  range  information  to 
clouds  and  fog.  A  lidar  inversion  was  conducted  without  prior  calibration  and  it  provided 
reasonable  cloud  structure  data. 

The  next  chapter  considers  the  integration  of  the  lidar  with  the  IRSTD. 
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V.  LIDAR  INTEGRATION  WITH  THE  NPS-IRSTD 


A .  CONCEPT 

The  general  concept  for  the  integration  of  the  lidar  with  the  NPS-IRSTD  system 
appears  to  be  fairly  simple  but  increases  in  complexity  with  a  detailed  analysis.  The  lidar 
data  would  be  used  initially  to  obtain  range  information  to  those  targets  identified  by  the 
NPS-IRSTD  but  a  secondary  aim  would  be  to  use  the  lidar  data  to  correct  the  IRSTD  image 
for  path  attenuation.  This  would  permit  enhancement  of  the  infrared  image  for  target 
identification  based  on  accurate  temperature  measurements.  The  primary  aim  of  the  thesis 
to  obtain  range  information  to  aunospheric  features,  such  as  clouds  and  fog,  using  the  lidar 
was  achieved,  but  in  the  process  enough  information  was  available  to  determine  the  actual 
visibility.  The  determination  of  the  visibility  could  be  very  important  for  defining  the  actual 
atmospheric  conditions  for  future  experiments.  The  ability  to  correct  the  NPS-IRSTD 
image  for  atmospheric  attenuation  was  not  an  easy  problem  as  there  was  no  direct 
correlation  available  between  visible  and  infrared  region  extinction.  A  possible  solution  has 
been  considered  which  would  require  additions  to  the  present  lidar  system. 

B.  LOCATION  CONSIDERATIONS 

A  number  of  details  will  need  to  be  considered  when  integrating  the  lidar  into  the 
existing  IRSTD  system;  physical  placement  of  the  two  systems  will  be  one  of  the  most 
important.  Firstly,  as  part  of  the  system  calibration  process,  the  ERSTD  and  the  lidar  will 
need  to  be  aligned  to  a  predetermined  range,  defined  by  operational  requirements  and  the 
IRSTD  system  characteristics.  The  placement  of  the  lidar  in  relation  to  the  IRSTD  will  be 
dependent  on  the  IRSTD’ s  location  on  board  the  ship  as  well  as  space  availability.  The 
following  are  a  number  of  placement  possibilities: 
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1 .  The  ideal  placement  of  the  lidar  would  be  inside  the  IRSTD  since  there  would  be  no 
misalignment  of  optical  axis  and  aiming  ambiguity.  It  would  not  be  possible  to  fit 
the  designed  lidar  into  the  IRSTD’s  enclosure  due  to  space  constraints.  It  may  be 
reconsidered  if  the  IRSTD  and  lidar  were  redesigned  to  be  integrated  in  the  same 
enclosure. 

2.  An  alternative  approach  would  be  to  have  the  lidar  placed  directly  above  or  below 
the  IRSTD.  This  would  maintain  the  rotational  mechanical  axis  alignment  at  the 
expense  of  the  vertical  pointing  accuracy.  In  the  limit,  as  the  range  becomes  much 
larger  than  the  separation  between  the  lidar  and  IRSTD,  the  vertical  pointing 
inaccuracy  becomes  negligible.  The  only  concern  would  be  the  supports  for  the 
two  sensors  and  whether  they  would  adversely  affect  each  other’s  field  of  view. 
This  would  require  an  operational  trade-off  study  to  determine  the  best  option. 

3.  The  most  probable  approach  would  be  the  placement  of  several  lidars  aboard  the 
ship.  Since  real  estate  on  the  ship's  superstructure  will  be  at  prime,  the  placement 
of  the  IRSTD  and  the  lidar  on  it  may  pose  some  problems  based  on  operational 
priorities.  Therefore  by  using  several  lidars,  each  lidar  could  cover  a  sector  of  the 
ship’s  perimeter.  The  coordination  of  the  lidars  with  the  IRSTD  may  be  complex 
and  the  pointing  accuracy  would  be  a  concern  which  would  require  detailed  study. 


There  is  one  location  consideration  which  only  applies  to  the  lidar  and  it  deals  with 
the  question  of  the  laser  eye  safety  requirements.  In  the  present  design,  laser  eye  safety 
was  designed  into  the  system  without  considering  an  operational  range  requirement  and  this 
permits  unrestricted  use.  The  range  requirement,  based  on  the  IRSTD  detection  range  and 
threat,  may  dictate  a  noneye  safe  laser  to  obtain  this  range  and  this  may  pose  operational 
problems  for  friendly  forces  and  restrict  its  use.  This  problem  is  similar  to  the  army’s 
concern  during  OPERATION  DESERT  STORM  with  the  increased  use  of  laser  range 
finders  on  the  battlefield  by  both  enemy  and  friendly  forces. 

The  IRSTD’s  maximum  detection  range  can  be  assumed  to  be  the  horizon  for  sea 
skimming  missiles  and  this  distance  can  be  expressed  in  terms  of  altitude  as: 


distance  (in  nautical  miles)  =1.14  s°-5 


(5-1) 
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where  s  is  the  altitude  in  feet  (Ref.  37].  If  the  IRSTD  is  placed  on  a  ship’s  superstructure 
at  60  feet  above  sea  level,  then  the  corresponding  range  is  16.4  kilometer  which  is  much 
larger  than  the  projected  range  of  the  designed  lidar.  Therefore,  to  achieve  this  range,  a 
noneye  safe  laser  will  need  to  be  employed. 

The  lidar’s  employment  may  be  restricted  to  those  ships  operating  on  the  perimeter 
of  a  carrier  group  and  further  restricted  from  use  in  designated  aircraft  inbound  and 
outbound  airspace. 

C.  RANGE  INTEGRATION 

The  determination  of  the  range  to  atmospheric  features  was  carried  out  by  using 
Klett’s  algorithm  as  outlined  in  Chapter  II.  The  time-dependent  data  obtained  from  the 
digital  oscilloscope  was  converted  into  range  data  using  the  speed  of  light.  The  lidar 
system  could  easily  be  converted  into  a  laser  range  finder  for  hard  targets  by  bypassing  the 
Klett  algorithm  and  using  a  reflected  signal  threshold  as  used  in  radar. 

Once  the  range  to  a  target  lias  been  determined,  the  display  of  this  information  on  the 
NPS-IRSTD  console  and  the  frequency  of  updating  need  to  be  addressed.  These 
considerations  are  dependent  on  the  type  of  target,  which  infers  its  speed,  the  IRSTD 
rotation  rate,  the  lidar’s  target  acquisition  time  and  human  visual  psychophysics. 

The  lidar  inversion  process,  besides  giving  range  information,  provides  enough 
information  to  determine  the  visibility.  The  determination  of  visibility  can  be  very  useful 
for  experimentation  purposes  as  it  provides  an  unambiguous  representation  of  the  weather 
conditions.  The  visibility  or  visual  range  can  be  determined  using  the  following  empirical 
formula  [Ref.  38]: 


Kv  X.  / 


(5-2) 
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where 


o(x)  is  the  mean  value  of  the  wavelength  dependent  extinction  coefficient  in 
kilometers  i, 

Rv  is  the  visibility  or  visual  range  in  kilometers. 

X  is  the  wavelength  in  nanometers,  and 

q  is  an  empirical  constant  which  is  defined  as  0.585  (Rv)1/3  for  Rv  <  6  kilometers  or 
=  1.3  for  average  seeing  conditions  [Ref.  39], 


It  has  been  noted  that  Equation  (5-2)  becomes  inaccurate  above  2  microns  producing 
a  significant  underestimation  of  the  extinction  coefficient  [Ref.  39],  Therefore,  visibility 
estimates  in  the  visible  spectrum  cannot  adequately  be  correlated  to  those  for  the  infrared 
spectrum. 


D.  INFRARED  EXTINCTION  COEFFICIENT 

Since  the  determination  of  the  infrared  extinction  becomes  difficult  when  using 
visible  spectrum  data,  the  best  approach  would  be  to  use  a  laser  which  operates  in  the 
region  of  the  infrared  sensor.  The  extinction  coefficient  can  be  directly  determined  using 
the  Klett  algorithm  and  the  infrared  image  can  be  directly  enhanced.  There  are  a  number  of 
disadvantages  in  this  approach  which  directly  influenced  the  choice  of  a  laser  for  this  thesis: 

1 .  The  power  law  backscattering  coefficient  and  extinction  coefficient  relationship  is 
not  as  strong  in  the  infrared  region  as  it  is  in  the  visible  spectrum  (as  discussed  in 
Chapter  II).  Therefore,  the  extinction  coefficient  is  harder  to  determine  in  the 
infrared  spectrum. 

2.  The  possibility  of  blinding  the  NPS-IRSTD  due  to  the  near-field  backscattercd 
radiation  is  a  concern  if  a  3-5  micron  laser  was  to  be  employed  in  the  lidar.  The 
use  of  a  laser  line  filter  on  the  IRSTD  may  be  considered  to  prevent  blinding  or 
damage. 

3 .  There  was  no  laser  available  at  the  school  operating  in  the  3-5  micron  region  for  the 
lidar  thesis.  Those  commercial  lasers  available  which  could  operate  in  this  region 
are  deuterium  fluoride,  helium  neon  (4-5  microns  reported  by  one  manufacturer), 
erbium: Y AG  (2.94  microns)  and  a  few  tunable  lead  salt  semiconductor  lasers  [Ref. 
40]. 
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1 .  Apparent  Temperature  Solution 

A  possible  solution  lor  the  determination  of  the  infrared  extinction  coefficient 
would  be  the  use  of  clouds  or  fog  infrared  temperatures.  In  the  infrared  spectrum,  satellite 
remote  sensing  uses  the  fact  that  clouds  can  be  approximated  as  blackbodies  with  a 
definable  temperature  and  this  approach  can  be  applied  in  this  application.  A  cloud’s 
temperature  can  be  obtained  from  the  NPS-IRSTD  but  this  temperature  would  be  the 
apparent  temperature,  since  there  would  be  path  attenuation  of  the  cloud  emission.  With  a 
modified  lidar,  the  path  temperature  up  to  and  including  the  cloud  or  fog  boundary  could  be 
determined.  The  mean  path  transmittance  can  then  be  calculated  assuming  that  the  cloud 
section  imaged  by  the  NPS-IRSTD  radiates  as  a  blackbody  into  a  hemisphere  with  the 
temperature  obtained  from  the  lidar. 

There  are  a  number  of  lidar  methods  used  to  measure  the  path  temperature 
profiles;  these  include:  [Ref.  41] 

1 .  Raman  scattering, 

2 .  differential  absorption,  and 

3 .  doppler  broadening. 

a .  Raman  Scattering 

Raman  scattering  is  the  weakest  of  the  scattering  phenomena  but  it 
provides  a  means  of  identifying  the  scattering  due  to  a  particular  molecule  regardless  of  the 
wavelength  of  the  incident  radiation,  which  makes  it  attractive  for  fixed  wavelength  lasers. 
Raman  scattering  is  inelastic  scattering  and  is  developed  from  the  quantum  theory  of 
scattering.  It  is  characterized  by  a  series  of  scattered  wavelengths  Xiaser,  Xiaser  ±  X,,  Xiaser 
±  X2,  etc.  where  Xj,  X2,  etc.  refer  to  the  energy  differences  in  the  molecular  species  in  the 
laser  beam  path.  The  longer  wavelength  shifted  lines  are  referred  to  as  Stokes  lines  while 
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the  shorter  wavelength  shifts  are  referred  to  as  anti-Stokes  lines.  Matched  pairs  of  shifted 
lines  (i.e.,  Xiaser  ±  X2)  have  roughly  the  same  intensity.  The  radiating  wavelength  will 
either  gain  or  lose  energy  from  the  incident  molecule’s  rotational  and/or  vibration  states. 
The  scattered  light  spectrum  will  have  a  strong  intensity  at  the  original  wavelength  and 
much  weaker  intensities  at  the  shifted  wavelengths.  This  molecular  dependent  shift  in  the 
wavelengths  is  called  the  characteristic  Raman  spectrum  and,  when  applied  to  remote 
sensing,  this  phenomenon  is  referred  to  as  Raman  scattering. 

The  pure  rotational  Raman  spectrum  has  been  used  for  lidar  applications 
because  the  magnitude  of  the  spectral  lines  is  two  orders  of  magnitude  greater  than  the 
vibrational-rotational  spectral  lines.  The  larger  spectral  magnitudes  make  the  pure  rotational 
Raman  spectrum  easier  to  measure  than  the  weaker  vibrational-rotational  spectrum  without 
loss  of  information.  The  intensity  of  the  pure  rotational  Raman  spectral  lines  can  be 
written  as  [Ref.  42]: 


I(J,T) 


lo-^gt 


BN 


kT 


2-(2J+  1)  S(J)  exp 


(5-3) 


where 


J  is  the  Stokes  and  anti-Stokes  Raman  rotational  quantum  number, 
Io  is  the  incident  intensity  of  the  laser  beam, 

Xj  is  the  Raman  spectral  wavelength  which  is  molecule  dependent, 
c  is  the  speed  of  light, 

gl  is  the  statistical  weight  due  to  nucleus  spin, 

B  is  the  rotational  constant  of  the  incident  molecule, 

N0  is  the  number  of  atoms  in  the  ground  state, 
k  is  Boltzmann’s  constant,  and 
T  is  the  temperature  of  the  molecule. 


The  product  of  the  S(J)  and  the  degeneracy  factor  (2J  +  1)  can  be  written 
for  the  Stokes  and  anti-Stokes  spectral  lines  as: 
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(Stokes  Branch)  (5-4) 


(2J  +  1 )  S(J) 


(J+  D(J  +  2) 
(2J  +  3) 


or 


(2J  +  l)S(J) 


J  (J  +  2) 
(2J-  1) 


(anti-Stokes  Branch).  (5-5) 


The  temperature  dependence  of  the  intensity  of  the  Raman  scattered 
wavelength  can  be  clearly  seen  from  Equation  (5-3).  There  are  a  number  of  Raman 
scattering  approaches  for  measuring  temperature  (Refs.  42,  43  and  44],  One  such 
approach  uses  two  rotational  Raman  scattering  molecular  lines  from  which  a  simple 
temperature  dependent  expression  can  be  written  (Ref.  42): 

R(D  =  exp(qr  +  P)  (5-6) 

where 

R(T)  is  the  ratio  of  the  Raman  scattered  intensities  or  received  powers, 
a  is  the  ratio  of  rotational  energies, 

T  is  the  temperature,  and 
P  is  ratio  of  the  difference  Raman  S(J)  functions. 

The  a  term  can  be  written  as: 


where 


_  ^rot( J 1 )  ~  Erotic) 
k 


(5-6) 


E  is  the  mean  rotational  energy  of  the  Raman  scattered  wavelength, 
J  is  the  Raman  quantum  number,  and 
k  is  Boltzmann’s  constant. 


The  p  term  can  be  expressed  as: 
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where  the  Raman  S(J)  function  can  he  determined  from  Equation  (5-4)  or  (5-5)  depending 
on  which  Raman  spectral  wavelength  is  observed. 

In  the  visible  spectrum,  the  popular  molecular  species  used  for  Raman 
scattering  are  oxygen  and  nitrogen  because  of  their  uniform  concentrations  in  the 
atmosphere.  The  only  problem  involved  in  this  approach  is  the  separation  of  the  Raman 
scattered  lines  without  significant  loss  of  signal  strength.  In  the  past,  this  has  been  a 
significant  problem  but  now  specially  developed  laser  line  notch  filters,  called  holographic 
notch  filters,  have  alleviated  this  problem.  These  filters  have  a  very  narrow  bandwidth 
which  make  them  ideal  for  the  removal  of  the  transmitting  laser  wavelength  from  the 
backscattered  signal.  An  optical  grating  could  then  be  used  to  observe  the  desired  Raman 
scattering  wavelength.  This  approach  to  temperature  measurement  would  be  easy  to 
implement  on  this  lidar  as  it  only  requires  modification  to  the  lidar  receiver. 
b .  Differential  Absorption 

The  differential  absorption  method  used  for  the  measurement  of 
temperature  profiles  requires  two  laser  wavelengths,  one  on  a  molecular  species  absorption 
wavelength  and  the  other  slightly  off.  The  ratio  of  the  backscattered  signals  is  used  to 
obtain  information  about  the  path  temperature  because  the  absorption  coefficient  is 
temperature  dependent.  The  two  wavelengths  used  need  to  be  close  enough  together  so 
that  the  transmittance  along  the  path  is  approximately  the  same  for  each  wavelength  and  the 
lasers  have  the  same  energy  output  to  remove  these  dependencies  from  the  ratio  of  the  two 
backscattered  signals.  The  typical  absorption  wavelengths  chosen  have  been  those  of 
water  vapor  or  oxygen  lines  (Ref.  45],  but  the  best  choice  would  probably  be  the  oxygen 
line  since  it  is  uniformly  mixed  in  the  atmosphere. 


In  order  to  use  this  approach  both  the  lidar  receiver  and  transmitter  would 
have  to  be  modified.  A  second  laser  with  the  same  energy  output  as  the  Nd.YAG  laser 
would  be  needed  as  well  as  a  modifications  to  the  receiver. 

c .  Doppler  Broadening 

The  doppler  broadening  approach  to  temperature  measurement  uses  the 
doppler  molecular  species  broadening  of  the  lidar  backscattered  frequency  spectrum.  This 
spectrum  is  composed  of  a  large  aerosol  component  and  weak  molecular  component.  The 
molecular  component  has  a  Gaussian  lineshape  which  is  related  to  temperature.  A  scanning 
Fabry-Perot  interferometer  could  he  used  to  measure  the  spectral  linewidth  between  the  half- 
width  half-maximum  points  and  thereby  obtain  a  temperature  (Ref.  46],  If  the  doppler 
broadening  is  known  for  a  reference  temperature,  then  the  following  expression  can  be 
written  [Ref  5]: 


AX(T)  —  AX^Tfgf) 


(5-8) 


where 

AX(T)  is  the  lidar  spectral  linewidth  at  temperature  T,  and 

AX(Tref)  is  the  lidar  spectral  linewidth  at  a  calibrated  reference  temperature  Tref. 

This  approach  would  be  difficult  to  implement  since  the  doppler 
broadening  varies  along  the  path  and  the  interferometer  may  not  be  able  to  scan  fast  enough 
to  determine  the  frequency  spectrum  along  the  path.  The  receiver  would  need  to  be  range 
gated  to  obtain  temperature  measurements  at  specified  ranges  such  as  a  cloud  boundary. 
The  boundary  temperature  would  be  used  for  the  apparent  temperature  solution  for  the  path 
average  infrared  extinction  coefficient  and  only  requires  additions  to  the  receiver.  This 
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approach  would  not  lend  itself  well  for  further  studies  of  path  temperatures  as  it  would  be 
very  slow  to  obtain  temperature  profiles  over  a  long  path. 

The  next  chapter  will  present  the  thesis  conclusions  and 
recommendations  for  future  further  lidar  work. 
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VI.  CONCLUSIONS  AND  RECOMMENDATIONS 


A  lidar  was  designed  and  manufactured  and  the  lidar  was  able  to  range  atmospheric 
features  as  demonstrated  in  Chapter  IV.  The  theoretical  maximum  range  for  the  lidar  is 
approximately  3.3  kilometers  by  night  and  0.9  kilometer  by  day.  These  ranges  are  well 
below  the  required  detection  range  for  integration  with  the  IRSTD,  which  was  determined 
to  be  16.4  kilometers  for  targets  on  the  horizon.  Never-the-less,  the  lidar  did  demonstrate  a 
ranging  capability  and  an  inversion  of  a  cloud  lidar  return  did  provide  an  estimate  of  a  cloud 
structure.  The  lidar  is  mounted  on  a  pedestal  which  allows  azimuth  and  elevation  pointing. 

There  are  many  improvements  to  be  made  on  the  lidar  which  were  discovered  during 
the  course  of  the  design,  manufacturing  and  evaluation  process: 

1 .  A  single  lens  was  used  to  expand  the  laser  beam  to  fill  the  telescope’s  secondary 
mirror  but  there  was  a  problem  with  overfill,  possibly  due  to  measurement  errors. 
A  double  lens  beam  expander  should  be  considered  to  provide  a  fine  adjustment  of 
the  beam  expander  focal  length  in  order  to  properly  fill  the  telescope’s  secondary 
mirror. 

2 .  The  laser  box  mirror  and  lens  mounts  used,  except  one,  were  all  non-precision 
mounts  and  to  get  the  correct  laser  beam  height  and  alignment  employed  a  trial  and 
error  method  for  the  placement  of  shims.  This  was  very  time  consuming  and,  if 
any  further  work  is  planned  for  the  laser  box,  adjustable  precision  optical  mounts 
should  be  purchased. 

3.  The  lidar  receiver’s  field  of  view  can  also  be  restricted  to  improve  the  daytime 
maximum  range  by  further  reducing  the  background  scene  noise.  This  would  not 
provide  any  increased  nighttime  range. 

4.  The  daytime  noise  levels  in  the  lidar  receiver  are  a  concern  and  need  to  be 
investigated  as  they  restrict  the  maximum  daytime  detection  range. 

5 .  If  the  lidar  daytime  receiver  noise  levels  can  be  reduced,  then  a  photomultiplier  tube 
amplifier  could  be  used  to  amplify  the  daytime  lidar  signals  to  further  increase 
detection  range. 
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6. 


A  logarithmic  amplifier  could  be  placed  at  the  output  of  the  lidar  receiver  to 
compress  the  lidar  signal.  This  amplifier  would  only  amplify  the  weak  long-range 
signals  while  leaving  the  near-field  return  in  its  original  form.  By  amplifying  the 
long  range  signals,  more  information  about  long-range  atmospheric  features  could 
be  obtained. 

7.  The  lidar  system  needs  to  be  calibrated,  as  discussed  in  Chapter  IV,  to  determine 
the  value  for  the  lidar  constant.  Without  a  calibration  of  the  lidar  system,  accurate 
path  extinction  coefficient  information  will  not  be  possible. 


This  lidar  is  viewed  as  a  stepping  stone  for  further  research  of  the  atmosphere.  A 
standard  lidar  inversion  technique,  the  “Klett”  method,  was  used  to  demonstrate  the  lidar’ s 
ability  to  provide  information  on  the  path  extinction  coefficient  both  using  theoretical  and 
real  lidar  data.  Further  work  and  refinement  is  needed  to  determine  the  best  lidar  inversion 
technique  and  then  to  apply  this  technique  to  the  lidar  measurements  taken  in  Chapter  IV  to 
see  if  better  results  can  be  achieved.  In  order  to  take  temperature  and  humidity  atmospheric 
measurements  in  the  future,  an  investigation  will  need  to  be  conducted  on  the  best  method 
to  be  used  from  those  described  in  Chapter  V  and  from  more  recent  sources.  If  the  Raman 
scattering  approach  is  to  be  used,  a  holographic  notch  filter  should  be  purchased. 

The  lidar  system  achieved  its  goal  to  provide  range  information  from  atmospheric 
features.  There  has  been  identified  a  number  of  areas  for  improvement  in  order  to  enhance 
the  lidar’s  accuracy  and  range  as  well  as  for  follow-on  research. 

This  thesis  has  been  a  challenging  and  enlightening  experience  in  the  area  of  lidar 
theory,  design  and  manufacturing. 
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APPENDIX  A 


ANALYSIS  OF  KLETT  EQUATION 


This  appendix  presents  an  analysis  of  the  Klett  equation  and  a  method  by  which  an 
expression  for  the  extinction  coefficient  boundary  value,  a(rm),  can  be  determined.  The 
Klett  equation,  as  discussed  in  Chapter  2,  produces  a  stable  solution  for  the  lidar  equation 
by  determining  <J(rm)  and  integrating  the  Klett  equation  backwards  from  rm  to  r  to  obtain 
the  extinction  coefficient  along  the  path,  o(r).  The  value  for  o(rm)  can  be  approximated  by 
relating  it  to  the  signal  strength  S(rin).  The  general  form  of  the  Klett  equation,  as  presented 
in  Chapter  2  as  Equation  (2-8),  is: 


CT(r)  = 


exp  [(S(r)  -  S(rm)yk] 


rim 

\  I  exp[(S(r')  -  S(rm)Vk]  dr’ 


(A-l) 


where 


a(r)  is  the  extinction  coefficient  along  a  path  at  range  r, 
a(r,n)  is  the  extinction  coefficient  boundary  value  at  range  rm, 

S(r)  is  the  logarithmic  range  adjusted  power  at  range  r,  S(r)  =  ln(  r2  P(r) ), 
S(rm)  is  the  logarithmic  range  adjusted  power  at  the  boundary  range  rm,  and 
k  is  the  exponent  used  in  the  relationship  (3(r)  =  Constant  a(r)C 


The  first  step  to  obtain  an  expression  for  a(r,n)  in  terms  of  S(rm)  is  to  make  a  variable 
substitution  into  Equation  (A-l)  and  then  solve  the  resulting  equation  as  a  linear  differential 
equation  as  shown  below: 
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1 .  A  new  variable  F(r)  and  its  derivative  are  defined  in  terms  of  the  denominator  of 
Equation  (A-l) . 


ffm 

F(r)  =  1  +  L^sl  exp  [(S(r')  -  S(rm))/k)  dr' 
d  F(r)  -  2  a(rm) 


dr 


exp  [(S(r)  -  S(rm))/k]  . 


(A-2) 

(A-3) 


2.  After  the  substitution  of  Equations  (A-2)  and  (A-3)  into  Equation  (A-l),  the  resulting 
expression  is: 


0|r)  . 

V  '  P(r\  <\r 


2  F(r)  dr 

3.  Equation  (A-4)  can  now  be  rewritten  as  a  linear  differential  equation. 

d  F(r)  ->  .  .  .  .. 

~7~~  +  T  a(r)  F(r>  =  0 
dr  k 


(A-4) 


(A-5) 


4.  The  solution  to  the  differential  equation  can  be  obtained  by  using  an  integrating 
factor.  The  choice  of  the  lower  limit  of  integration  for  the  integrating  factor  is 
arbitrary. 


Integrating  Factor:  exp  |  2-  |  o(r)  dr 


f 


(A-6) 


JL 

dr 


F(r)  exp  |  2-  j  <J(r)  dr  | 


=  0 


*rm 

r  ~\\ 

1 

,  r' 

^  F(r)  exp 

kj  °(r)drJj 

J  ro 

dr’  =  0 


(A-7) 


(A-8) 


F(rm)  exp  f- 


L  rro 

kjr,  °(r)dr 

-  F(r0)exp  JJ  a(r)dr 

=  0 


(A-9) 
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5.  The  selection  of  the  limit  of  integration,  r',  should  now  be  considered  in  order  to 
simplify  the  analysis.  By  selecting  r'  =  rm  and  noting  that  F(rm)  =  1  as  determined 
from  Equation  ( A-2),  the  resulting  expression  for  Equation  (A-9)  becomes: 


exp  2.  a(r)  dr 

•^m 


1 

F(r0) 


(A- 10) 


6.  If  the  natural  logarithm  is  taken  of  Equation  (A- 10)  and  the  limits  of  integration  for 
the  integral  are  switched,  the  resulting  expression  is: 

a(r)dr  =  ^  ln(F(r0))  (A-ll) 

r0 


By  substituting  Equation  (A-2)  for  the  expression  for  F(ru)  into  Equation  (A-l  1),  the 
resulting  expression  is  the  same  as  that  used  by  Klett  [Ref.  20): 


rm 

I  a(r)  dr’  =  ^  In 

1  +  exp[(S(r')-S(rm))/k]dr' 

Jto 

(A- 12) 


or,  using  a  more  compact  notation: 


o(r)  dr'  =  ~  In  ( 1 


(A- 13) 


where 


I  =  (rm-  r0)  L 


exp[(S(r’)  -  S(rm))/k]  dr' 


(A- 14) 
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and 


2  o( rm)  (rm  -  r„) 


(A-15) 


Using  the  assumption  that  the  average  extinction  coefficient  over  the  range  0  to  rm  is 
equal  to  that  over  the  range  r()  to  r,n,  which  is  reasonable  if  rG  «  rm,  additional  progress 
can  be  made  on  Equation  (A- 13). 

1 .  The  average  path  extinction  coefficient  is  approximated  as: 


Jfrn>  f1  n 

a(r)dr  =  — ! — - 

o  (Fm  * r<>)  JrQ 


a(r)  dr 


( A- 16) 


2.  The  expression  for  S(rm)  is  derived  from  Equations  (2-1)  and  (2-3). 


/ 


S(rm)  =  In 


P0  ^  A  Constant  a(rm)k  exp  -  2 1  a(r)  dr 


frm 

JO 


Jprm 

a(r) 

o 


V 


dr 


(A- 17) 


(A- 18) 


S(rm)  =  C  +  k  In  (a(rm))  -  2  o(r)  dr 


Jrrm 

o 


(A- 19) 


3.  The  integral  term  in  Equation  (A- 19)  can  be  rewritten  using  the  equivalence 
relationship  from  Equation  (A- 16). 


S(rm)  =  C  +  k  In  (o(r,n))  - 


tTm 

2rm 

(tm  '  r0)  I 

^rO 


a(r)  dr 


(A-20) 
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4.  The  integral  term  in  Equation  ( A-20)  can  be  removed  by  the  substitution  of  Equation 
(A- 13)  and  the  resulting  expression  is: 


S(rm)  -  C  =  k  In  (<T(rin))  -  k.  ln  ( 1  +  I  flj 

O in  '  ro)  4 


(A-21) 


5.  By  adding  to  both  sides. 


ln(  — r*"  r°) 


(A-22) 


the  linal  expression  for  o(rm)  in  terms  of  S(rm)  is  now  in  the  form  presented  as 
liquation  (2-9). 


S(rm) '  C  f  In  (2  '  r“>)  =  In  (On)  ■  j In  ( 1  +  I  £lm) 


(A-23) 
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APPENDIX  B 


THEORETICAL  LIDAR  INVERSION 

A.  INTRODUCTION 

This  appendix  examines  the  theoretical  aspects  of  the  “Klett”  method  for  lidar 
inversions.  For  the  analysis,  eight  known  extinction  profiles  (o(r))  are  used  to  generate  a 
S(r)  or  logarithmic  range  adjusted  power  profile.  Then  using  the  Klett  inversion  method, 
the  S(r)  profile  was  inverted  to  obtain  the  extinction  coefficient  profile.  A  comparison  of 
the  inverted  to  the  known  extinction  profile  was  made  using  low  and  high  range  resolution 
data.  A  brief  second  analysis  was  conducted  on  selected  extinction  profiles  by  adding 
Gaussian  white  noise  on  top  of  the  received  lidar  signal  to  determine  the  effect  of  noise  on 
the  inversion  process. 

B.  EXTINCTION  PROFILES 

In  order  to  test  the  Klett  inversion  algorithm,  a  number  of  extinction  profiles  were 
obtained  fiom  Mr.  William  Lent/,,  a  physicist  at  the  Naval  Postgraduate  School,  who  has 
conducted  previous  work  in  this  field  and  has  worked  with  Mr.  Klett.  The  extinction 
profiles  chosen  were  used  by  Mr.  Lentz  in  his  work  in  lidar  inversions  and  are  typical  of 
extinction  profiles  which  could  be  encountered. 

The  first  set  of  extinction  profiles,  curves  a  and  b  shown  in  Figure  14,  deal  with  the 
simplest  case,  in  which  the  extinction  coefficient  is  constant.  The  next  set  of  extinction 
coefficients,  curves  c  and  d  shown  in  Figure  15,  shows  a  lidar  profile  going  from  a  low  to 
high  visibility  (e.g.,  out  of  fog  bank  or  cloud).  The  constant  extinction  coefficient  out  to 
255  meters  simulates  the  Klett  assumption  that  the  extinction  coefficient  from  range  0  to  r0 
is  constant  and  approximately  equal  to  the  extinction  coefficient  at  range  r0,  as  described  in 
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Appendix  A.  The  range  r0  is  the  minimum  range  at  which  the  lidar  transmitter  and  receiver 
fields  of  view  overlap  and  meaningful  lidar  data  can  be  taken. 


The  next  set  of  extinction  profiles,  curves  e  and  /  shown  in  Figure  16,  are  the 
inverse  of  Figure  15.  In  this  case,  the  lidar  profile  simulates  going  into  a  cloud  or  fog 
while  the  lidar  signal  fades  into  noise.  As  in  the  Figure  15,  the  Klett  assumption 
concerning  the  constant  nature  of  the  extinction  coefficient  over  the  range  0  to  r0  is  used. 
The  final  set  of  extinction  profiles,  curves  g  and  h  shown  in  Figure  17,  consider  the  case 
where  the  lidar  is  able  to  penetrate  through  a  cloud  or  fog  bank  and  obtain  data  from  the 
other  side. 


Range  (m) 

Figure  14.  Constant  Extinction  Coefficient  curves  a  and  6. 
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Extinction  Coefficient  (kmA- 1 )  frs  Extinction  Coefficient  (kmA- 1 ) 
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ure  15.  Exponentially  Decreasing  Extinction  Coefficient  curves  c  and  d . 


Figure  16.  Exponentially  Increasing  Extinction  Coefficient  curves  e  and /. 


3 


Range  nn) 


Figure  17.  Dynamic  Range  Extinction  Coefficient  curves  g  and  h  . 


C.  ANALYSIS 

The  lidar  theoretical  inversion  program,  listed  in  Section  D,  was  written  in  the 
MATLAB  numeric  computation  software  programming  language.  The  process  used  to 
invert  a  lidar  return  was: 

1 .  Convert  the  known  extinction  coefficient  profile  into  a  S(r)  profile  using  the  general 
form  of  Equation  (A- 17)  rewritten  as: 


S(r)  =  In  jp0  &  A  Constant  j  +  k  In  (a(r))  -  2 1  a(r)  dr  •  (B-l) 

Jo 

The  ln(P0  ct/2  A  Constant)  term  is  a  constant  and  can  be  set  to  any  arbitrary  value. 
The  determination  of  this  value  is  unimportant  since  it  will  be  canceled  out  during  the 
inversion  process. 

2.  Invert  the  generated  S(r)  profile  using  the  selection  criteria  for  the  extinction 
coefficient  boundary  value  (a(rm))  described  in  Chapter  II. 


74 


3.  Compare  the  known  extinction  coefficient  profile  to  inverted  extinction  coefficient 
profile. 

The  lidar  inversion  algorithm  was  applied  to  simulated  low-resolution  (7.5  meter) 
and  high-resolution  (0.15  meter)  range  data.  The  range  resolution  is  determined  by  the 
sampling  rate  used  to  convert  the  analog  backscattered  return  to  digital  form.  The  low- 
resolution  case  was  provided  by  Mr.  William  Lentz  with  the  sample  extinction  profiles 
where  the  sampling  rate  used  was  20  megaHertz  [Ref.  9].  The  high-resolution  case  was 
analyzed  because  of  its  important.,  m  the  actual  lidar  inversion  to  be  theoretically  achieved 
in  this  thesis.  The  backscattered  return  was  initially  planned  to  be  digitized  using  a 
digitizing  oscilloscope  which  had  a  sampling  rate  of  1  gigasample  per  second.  The  time 
domain  sampling  interval  is  1  nanosecond  between  samples  which  converts  into  0.3  meters 
between  samples,  but  this  is  round-trip  time  and  it  should  be  divided  by  two  to  get  the 
range  from  which  the  backscattered  signal  originated.  Therefore,  the  actual  range 
resolution  is  0.15  meters. 

The  lidar  inversions  of  the  S(r)  files,  generated  from  the  extinction  coefficient  curves 
a-h ,  shown  at  Figures  15-17,  were  conducted  using  both  the  low  and  high-resolutions.  It 
was  rioted  during  the  lidar  inversion  analysis  that  most  of  the  inverted  profiles  were  very 
close  to  the  actual  extinction  profiles.  When  initially  displaying  the  data,  it  appeared  that 
one  or  more  profiles  are  missing  but  in  fact  they  were  actually  coincident.  Therefore,  it 
was  decided  to  display  only  the  differences  between  the  inverted  and  actual  extinction 
coefficient  profiles  so  as  not  to  confuse  the  reader.  The  following  observations  were  made 
from  each  lidar  inversion: 
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Extinction  Deviation  (kmA-l ) 


1.  Figure  18.  This  figure  shows  the  difference  between  the  low  and  high-resolution 
inversions  of  the  S(r)  file,  generated  from  curve  a  from  Figure  14.  and  the  actual 
curve  a..  The  high  range  resolution  inversion  followed  the  actual  extinction  profile 
much  better  than  the  low-resolution  case.  The  low-resolution  inversion  deviation  is 
mainly  due  the  numerical  integration  methods  which  approach  the  continuous  integral 
solution  when  the  interval  width  becomes  small  enough. 
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Figure  18.  Extinction  Coefficient  deviations  from  the  actual  values  (Figure  14  curve  a  )  by  the 
the  high  and  low-recolution  inversions. 
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Extinction  Deviation  (kmA-l) 


2.  Figure  19.  This  figure  shows  the  difference  between  the  low  and  high- resolution 
inversions  of  the  S(r)  file,  generated  from  curve  b  of  Figure  14,  and  the  actual  curve 
b.  The  high  range  resolution  inversion  followed  the  actual  extinction  piofile  better 
but  there  was  no  real  advantage  for  using  the  high-resolution  to  invert  the  S(r)  file 
when  the  low-resolution  inversion  worked  just  as  well.  Therefore  for  high  constant 
visibility  conditions,  low- resolution  data  acquisition  will  provide  accurate  results. 
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Figure  19.  Extinction  Coefficient  deviations  from  the  actual  values  (Figure  14  curve  b )  by  the 
the  high  and  low-resolution  inversions. 


Extinction  Deviation  (kmA-l ) 


3.  Figure  20.  This  figure  shows  difference  between  the  low  and  high-resolution 
inversions  of  the  S(r)  file,  generated  from  curve  c  from  Figure  15.  and  the  actual 
curve  c.  As  noted  in  Figure  IX.  the  high-resolution  inversion  provided  a  more 
accurate  representation  of  the  actual  Extinction  Coefficient  profile.  It  was  noted  that 
the  maximum  extinction  deviation  is  much  larger  than  in  the  two  previous  figures. 
Figure  18  and  19. 


range  (m) 

Figure  20.  Extinction  Coefficient  deviations  from  the  actual  values  (Figure  15  curve  c )  by  the 
the  high  and  low-resolution  inversions. 


78 


Extinction  Deviation  (kmA-l ) 


4.  Figure  21.  This  figure  shows  me  difference  between  the  low  and  high-resolution 
inversions  of  the  S(r)  file,  generated  from  curve  d  from  Figure  15,  and  the  actual 
curve  d.  As  noted  in  Figure  19,  there  was  no  advantage  in  using  the  high-resolution 
inversion  over  the  low-resolution  inversion,  since  the  maximum  deviation  was  in  the 
order  of  1()-6  km-i. 


range  (m) 

Figure  21.  Extinction  Coefficient  deviations  from  the  actual  values  (Figure  15  curve  d )  by  the 
the  high  and  low-resolution  inversions. 
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Extinction  Deviation  (kmA-l) 


5.  Figure  22.  This  figure  shows  the  difference  between  the  low  and  high-resolution 
inversions  of  the  S(r)  file,  generated  from  curve  e  from  Figure  16,  and  the  actual 
curve  e .  The  large  deviation  was  noted  in  both  the  high  and  low-resolution 
inversions.  This  is  due  to  the  Low  Visibility  algorithm.  The  algorithm  becomes 
sometimes  unreliable,  as  also  noted  by  Klett.  The  simple  assumption  that  the  mean 
of  the  extinction  coefficient  is  equal  to  the  boundary  value  extinction  coefficient. 
Equation  (2-16),  is  always  not  valid  when  there  is  a  large  dynamic  range  of  the 
extinction  coefficient. 


range  (m) 

Figure  22.  Extinction  Coefficient  deviations  from  the  actual  values  (Figure  16  curve  e  )  by  the 
the  high  and  low-resolution  inversions. 
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6.  Figure  23.  This  figure  shows  the  difference  between  the  low  and  high-resolution 
inversions  of  the  S(r)  file,  generated  from  curve/  from  Figure  16,  and  the  actual 
curve/.  The  large  deviation  noted  in  Figure  22  was  not  observed  in  Figure  23  which 
was  consistent  since  there  was  a  small  dynamic  range  for  the  Extinction  Coefficient. 


range  (in) 

Figure  23.  Extinction  Coefficient  deviations  from  the  actual  values  (Figure  16  curve /)  by  the 
the  high  and  low-resolution  inversions. 
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7.  Figure  24.  This  figure  shows  difference  between  the  low  and  high-resolution 
inversions  of  the  S(r)  tile,  generated  from  curve  g  from  Figure  17,  and  the  actual 
curve  g.  The  large  deviation  at  700  meter  range  previously  described  in  Figure  22 
for  the  behavior  of  the  Low'  Visibility  algorithm  was  again  observed  here.  The  effect 
of  increasing  the  range  resolution  is  to  increase  slightly  the  range  from  which 
accurate  Extinction  Coefficient  data  can  be  obtained. 


range  (m) 

Figure  23.  Extinction  Coefficient  deviations  from  the  actual  values  (Figure  17  curve  g  )  by  the 
the  high  and  low-resolution  inversions. 
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Extinction  Deviation  (kmA-l) 


8.  Figure  25.  This  figure  shows  the  difference  between  the  low  and  high-resolution 
inversions  of  the  S(r)  file,  generated  from  curve  h  from  Figure  17,  and  the  actual 
curve  h.  The  same  lack  of  effect  by  the  Low  Visibility  algorithm  was  noted.  The 
maximum  deviation  is  in  the  order  of  KM. 
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Figure  25.  Extinction  Coefficient  deviations  from  the  actual  values  (Figure  17  curve  h  )  by  the 
the  high  and  low-resolution  inversions. 


It  was  not  the  aim  of  this  thesis  to  improve  the  Klett  algorithm.  However,  it  may  be 
noted  that  the  unreliability  of  the  Low  Visibility  algorithm  may  not  be  as  much  of  a  problem 
as  first  thought.  The  main  problem  occurs  near  the  boundary  range  where  the 
backscattered  return  will  be  very  close  to  the  noise  floor.  Therefore,  the  last  part  of  the 
lidar  inversion  may  not  provide  any  useful  data  due  to  noise  corruption  and  could  be 
discarded. 
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The  final  part  of  the  lidar  inversion  analysis  involves  the  addition  of  noise  to  the 
received  backscattered  signal.  The  photon  noise  statistics  are  a  poisson  process  but,  in  this 
case,  since  the  lidar  receiver  observation  time  is  long,  the  Central  Limit  Theorem  can  be 
applied.  The  result  is  that  the  background  noise  can  be  approximated  by  a  Gaussian 
process  which  can  be  easily  implemented.  The  total  signal  at  the  output  of  the  receiver  is  a 
summation  of  the  actual  lidar  backscattered  signal  P(r)  and  the  background  noise  power. 
An  additive  noise  magnitude  was  arbitrarily  set  to  20%  of  the  value  for  P(rm)  in  the  absence 
of  actual  data. 

A  constant  (curve  a ,  Figure  14),  dynamic  range  (curve  h  ,  Figure  17)  and 
exponentially  increasing  (curve  Figure  16)  extinction  coefficient  profiles  were  used  in  the 
noise  analysis.  The  exponentially  decreasing  extinction  coefficient  profile  was  found  to 
have  the  same  noise  characteristics  as  the  exponentially  increasing  case  and  was  not  used. 
The  low-resolution  profiles  were  used  for  the  analysis  because  the  effects  were  easier  to 
observe. 

The  noisy  inverted  extinction  coefficient  profiles  generated,  from  curves  a,f  and  h, 
arc  in  Figures  25,  26  and  27,  respectively.  There  is  an  effect  of  the  simulated  background 
noise  on  the  inversion  but  it  influences  some  profiles  more  than  others.  In  general,  the 
suppression  of  the  background  noise  will  be  important  in  order  to  reach  the  theoretically 
predicted  maximum  ranges  of  the  lidar  in  Chapter  II.  There  are  a  number  of  approaches, 
such  as  frequency  spectrum  analysis  based  on  statistical  signal  processing  theory  and  time 
domain  analysis,  which  could  be  employed.  The  time  analysis  could  easily  be  employed 
by  averaging  a  number  of  lidar  returns  over  a  short  period  of  time  using  the  sample  and 
average  capability  on  the  digitizing  oscilloscope.  This  assumes  that  the  atmosphere  is 
relatively  calm  and  the  extinction  coefficient  profile  is  constant  over  a  short  period  of  time. 
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igure  25.  Noisy  Constant  Extinction  Coefficient  Inversion  using 
Figure  14  curve  a  profile. 
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Figure  26.  Noisy  Exponentially  Increasing  Extinction  Coefficient 
Inversion  using  Figure  16  curve /  profile. 
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Figure  27.  Noisy  Dynamic  Range  Extinction  Coefficient  Inversion 
using  Figure  17  curve  h  profile. 


D.  MATLAB  LIDAR  INVERSION  PROGRAM 

This  section  contains  the  MATLAB  code  used  to  perform  the  lidar  inversions. 

%  Theoretical  Lidar  Inversion  of  known  extinction  coefficient  profiles 
%  using  Klett's  Technique. 

%  This  program  is  robust  enough  to  use  high  or  low  resolution  (time/range  data) 

%  Author  -  Major  M.M.  Regush 
%  Date  -  26  May  1993 


clg 

clear 


% 

%  Load  Known  extinction  coefficient  files  low  or  high  range  resolution 
%  Constant,  Exponentially  Increasing/Decreasing  or  Dynamic 
%  range  extinction  profiles 
% 
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name=input( 'Input  Data  file  name  -  Vs'); 
eval(['load  '.name)) 
sigma=eval([namej); 
sigma=sigma'; 

%  Note  extinction  coefficient  file  starts  at  range  =  rQ. 

%  Therefore  must  add  one  more  extinction  value  at  the  start  to  obtain 
%  the  range  r  =  0  to  750  meters 
sigma=[sigma(l)  sigma); 

(p,q)=size(sigma); 


%  Definitions 

k=l;  %  power  law  relationship  exponential  between  the 

%  backscattered  &  extinction  coefficients 
range=linspace(0,7 50,q)/  l(KK); 

[m,n]=size(range); 

h=abs(range(2)-rangc(l));  Range  resolution 


% 

%  Determination  of  S(r)  -  Log  Adjusted  Power  P(r) 

% 

0  =  1;  %  Arbitrary  Constant  from  equation  (B-l) 

%  Integration  of  Extinction  Coefficient  using  trapezoidal  rule 

total=0; 

for  x=l:n 

total=total+sigma(x ) ; 
sint(x)=total-.5*(sigma(x)+sigma(  1 )); 
end 

sint=sint*h; 

Sr=Cl+k*log(sigma)-2*sint; 

%  Remove  the  range  =  0  m  point  as  it  is  no  longer  needed 
Sr=Sr(2:n); 

So=Sr(l); 
range=range(2:n); 
ro=range(l); 
sigma=sigma(2:n); 
n=n- 1 ; 


% 

%  Noise  added  to  the  lidar  backscattered  return,  if  desired. 

% 

noise=input('Add  simulated  Gaussian  background  noise  (y/n)  ->  Vs'); 
if  noise  =  'y' 

%  Make  P(r)  file 

Pr=exp(Sr)7(range.A2); 

rand(’normal')  %  Gaussian  distribution 

nr=rand(l,q-l); 

dispC  Arbitrarily  set  max  noise  magnitude  to  20%  of  P(rm)  magnitude') 
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nr=nr.*(Pr(q- 1  )*.2)/max(m ); 

Pr=Pr+nr; 

Sr=log(Pr.*(range.A2));  7<  New  modified  S(r)  file 

end 


% 

%  Lidar  Signal  Inversion  Analysis  -  sigma(r) 

% 

Sm=Sr(n); 
rm=max(  range); 

7c  Integration  of  exp(Sr-Sm ) 

snevv=exp((Sr-Sm)/k); 
total=(); 
for  x  =  l:n 
y=n-x+l; 

total=lotal+snew(y); 

Sint(y)=total-.5'l‘(snew(y)+snew<n)); 

end 

Sint=Sinl+h; 

%  Boundary  Value  Sigma  Determination 
flag=();  7c  Used  for  sigma(rm )  algorithms 
I=Sint(l)/(rm-ro); 

Gm=(Sm-Cl)/k  +  log(2*(rm-ro)/k); 
dispC  ') 

dispf'Use  High  Vis  Algorithm  to  Determine  sigma(rm)’) 
dispC  ') 

sro=0;  %  Initial  Condition  for  sigma(ro) 

Gmp=Gm  +  2*ro*sro/k; 

7c  Initial  Check  of  High  Vis  Algorithm 
if  exp(-Gmp)  <=  (I+.01) 
flag=  I ; 
end 

dsro=abs(Gm); 
while  dsro  >  ,0(KK)1 
Grnp=Gm  +  2*ro*sro/k; 
srm=k*.5/((rm-ro)*(exp(-Gmp)-I)); 
srol=snew(l)*srm/(l+2*srm/k*Sint(  1)); 
dsro=abs(sro-srol ); 
sro=sro  1 ; 
end 


7c  Final  Check  Of  High  Vis  Algorithm 
ifexp(-Gmp)  <=  (I+.01)  I  srm  <=  .01  I  (sro/srm)  >=  50 
fiag=l; 
end 
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if  flag  ==  1  &  I  >  1 

dispC***  High  Vis  Algorithm  Failed  Use  Low  Vis  Algorithm  ***') 
dispC  ') 

srm  1  =  1;  %  Initial  Condition  for  sigmat  rm ) 

gM=2*srm  1  *(rm-ro)/k; 
dgM=gM;; 
while  dgM  >  .1)01 
gMl=log(l  +  gM*I); 
dgM=abs(gM-gM  1 ); 
gM=gMl; 
end 

srm=gM*k*.5/(rm-ro);  %  Boundary  Value  sigma(rm) 
elseif  flag  ==  1  &  I  <  1 

dispC***  Low  Visibility  Algorithm  Failed  Use  Default  Algorithm  ***') 
dispC  ') 

gMc=(rm-ro)/(ro*I); 

srm=gMc*k*.5/(rm-ro);  %  Boundary  Value  sigma(rm) 
end 


% 

%  Lidar  Inversion  Eqn 

% 

snewp=(srmA(- 1  )+2/k*Sint); 
sn=snew.*snewp.A(- 1); 


% 

%  Plots  of  S(r)  file  and  actual  &  inverted  extinction  profile  results 
% 

subplot(21  l),plot((range-ro)*  l(KK),(Sr-So)), 

title('Log  Range- Adjusted  Power  from  Extinction  Coefficient  Profile’), 
xlabel('r-ro  (my),ylabel('S(r)  -  S(ro)') 
grid 

subplot(212),plot((range-ro)*l(XX),sn,,--',(range-ro)*100(),sigma,’-'), 
title('Actual  (  )  &  Inverted  Extinction  Profiles  (--)  Comparison'), 

xlabel(’r-ro  (m)'),ylabel('Extinction  (kmA-l)') 
grid 
pause 
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% 

%  Numerical  Comparisons 
% 

dispONumerical  Comparisons  of  Actual  &  Inverted  Extinction  Profiles'); 
fprintfOInverted  Extinction  Coefficient  (s>  rm  %g  kmA-l\ri,srm) 
fprintf ('Known  Extinction  Coefficient  <a>  rm  %g  kmA-l\n',sigma(n)) 
sbve=(srm-sigma(n))/sigma(n)+  1(K); 

fprintf('Extinction  Coefficient  Boundary  Value  Error  %g  Percent  \n',sbve) 
fprintf('Average  Inverted  Extinction  Coefficient  %g  kmA-l\n',mean(sn)) 
fprintf('Average  known  Extinction  Coefficient  %g  kmA-l\n',  mean(sigma)) 
ase=(mean(sn)-mean(sigma))/mean(sigma)*  1 00; 
fprintf('Average  Extinction  Coefficient  Error  %g  Percent  \n’,ase) 
fprintf('Delta  of  Extinction  Coefficients  @  ro  %g  \ri,  abs(sigma-sn(  1 ))) 
fprintf('Delta  of  Extinction  Coefficients  @  rm  %g  \n',abs(sigma(n)-sn(n))) 
mdbkis=max(abs(sigma-sn)); 

fprintfCMax  Difference  between  Extinction  Coefficients  @  %g  \n',mdbkis) 

%  Visual  Range  Determination  using  empirical  formula  (see  Chapter  V) 
disp('  ’) 

disp('Determination  of  Visibility  ) 
avis=3.91/mean(sigma)*(550/532)A1.3; 
ivis=3.91/mean(sn)*(550/532)A1.3, 
fprintf  ('Actual  Visibility  =  %g  km  \n’,avis); 

fprintf('Visibility  determined  from  Inversion  Process  =  %g  km  \n',ivis); 
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